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ABSTRACT 

We used the Plateau de Bure Interferometer (PdBI) to map a sample of 15 submillimetre galaxies (SMGs) in the COSMOS field at the 
wavelength of 1.3 mm. The target SMGs were originally discovered in the James Clerk Maxwell Telescope (JCMT)/AzTEC 1.1 mm 
continuum survey at S/N Umm - 4 - 4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations 
of these sources. The angular resolution of our observations, ~ 178, allowed us to accurately determine the positions of the target 
SMGs. Using a detection threshold of S/N L3mm > 4.5 regardless of multiwavelength counterpart association, and 4 < S/N L3mm < 4.5 
if a multiwavelength counterpart within 175 is also present, the total number of detections in our survey is 22. The most significant 
PdBI detection of S/Nx 3mm = 10.3 is towards AzTEC 19. Three of our detected SMGs (AzTEC21, 27, and 28; which corresponds to 
20%) are marginally resolved at our angular resolution, and these sources are found to have elongated or clumpy morphologies and/or 
multiple components. Using optical to near-infrared photometric redshifts, available spectroscopic redshifts, and redshifts estimated 
from the radio-to-submm spectral index we infer a median redshift of z = 3.20 + 0.25 for our sample. To study the overall multiplicity 
and redshift distribution of flux-limited samples of SMGs we combined these sources with the 15 brightest JCMT/AzTEC SMGs 
detected at 1.1 mm, AzTEC 1-15, and studied previously. This constitutes a complete, flux- and S/N-limited 1.1-mm selected sample. 
We find that the median redshift for the 15 brightest JCMT/AzTEC SMGs (z - 3.05 + 0.44) is consistent with that for AzTEC16-30. 
This conforms to recent observational findings that SMGs do not exhibit any significant trend between the redshift and (sub)mm flux 
density. For the combined AzTECl-30 sample we derive a median redshift of z = 3.17 + 0.27, consistent with previous results based 
on mm-selected samples. We further infer that within the combined AzTEC 1-30 sample ~ 25 + 9% of sources separate into multiple 
components. 
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1. Introduction 

When the first extragalactic submillimetre continuum surveys 
were carried out at the end of the 1990s, a population of 
heavily dust-obscured (i.e. optically faint) galaxies at high 


* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany), and IGN (Spain). 


redshift was disco vered (ISmail et al. 19971 |Hughes et al. 1 998; 
Barger et al. 1998). These sources are generally referred to 
as submillimetre galaxies or SMGs (see IBlain et al. 20021 
|Casey et al. 2014| for reviews). 

The bulk of SMGs are observed at redshifts z - 2 - 3 (e.g. 
Chapman et al. 20 05 War dlow et al. 20111 ILindner et al. 20 lH 
Casey et al. 2013| |Simpson et al. 2014| lUmehata et al. 20141) . 
However, the number of known high-redshift (z > 4) 
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SMGs has increased considerably in the past few years (e.g. 


Schinnerer et al. 2008; Daddi et al. 2009a,b; Coppin et al. 2009, 

Riechers et al. 2010; Capak et al. 2011; Smolcic et al. 2011, 


Eox et al. 2011; Combes et al. 2012; Walter et al. 2012. 

Swinbank et al. 2012; WeiB et al. 2013; Smolcic et al. 2015). 
The SMG with the highest spectroscopically confirmed redshift 
currently known is HFFS3 at z = 6.34 {Riechers et al. 2013), 
which demonstrates that these sources were already present 
when the universe was only < 890 Myr old (see the end of this 
section for our adopted cosmology). 

Submillimetre galaxies have parent dark matter haloes - 
i.e. the sites of galaxy formation originating in the dark-matter- 
dominated density perturbations in the early universe (e.g. 
Benson 2010) - with characteristic masses of ~ 10 12 - 10 13 M© 

|Blain et al. 2004; Swinbank et al. 2008; Hickox et al. 2012). 
The physical properties of SMGs are found to be extreme. 
In particular, their very high infrared (IR; 8-1 000 fi m) lumi¬ 
nosities of L ir ~ 10 12 - 10 13 F 0 are indicative of extreme 
star formation rates (SFRs) of ~ 100 - 1000 M© yr _1 , 

making SMGs the most intense known starbursts in the uni¬ 
verse. Observations of CO rotational transitions with upper 
rotational-energy levels of J u = 2 - 7 suggest H 2 gas masses 

of M h , ~ 10 10 - 10 11 M 0 in SMGs (e.g. [Greve et al. 2005} 
Tacconi et al. 2006; Bothwell et al. 2013 ), while CO(/ =1-0) 
observations yield gas masses up to several times 10 11 M© 
([ivison et al. 2011; Riechers et al. 2011). Submillimetre galax- 
ies are therefore among the most gas-rich systems in the uni¬ 
verse. For instance, the median Mh 2 value of 3.0 ± 1.6 x 10 10 M© 
(within a ~ 2 kpc radius) derived for SMGs studied by 

Greve et al. (2005) is about four times higher than in the 
most luminous local ultraluminous IR galaxies or UFIRGs 
(Solomon et al. 1997). These authors also estimated that the 
typical gas-consumption timescale in SMGs is > 40 Myr, 
but they noted that if feedback processes slow down the star 
formation activity (i.e. negative feedback such as radiation 
pressure acting on dust, stellar winds, outflows, supernovae, and 
the associated turbulence), the above timescale can be signifi¬ 
cantly longer. The derived stellar masses in SMGs are typically 

in the range M* ~ 10 11 - 10 12 M© (e.g. [Borys et al. 2005 

Dye et al. 2008; Wardlow et al. 2011; Hainline et al. 2011, 


Michalowski et al. 2012; Simpson et al. 2014). While some 

authors suggest that SMGs might predominantly constitute the 
high-mass end of the star-forming galaxies’ main sequence (the 
M+- SFR relationship) at z > 2 (Michalowski et al. 2012), a fair 
fraction certainlv lies above it (e.g. Daddi et al. 2009a|. 

Since SMGs are found to have very high SFRs, the 
question then arises as to which physical process(es) are 
responsible for these rates. It has been suggested that galaxy 
mergers can trigger a significant burst of star formation (e.g. 
Barnes & Hernquist 1991). The basic idea behind this is that, 
when dynamical friction within a parent halo causes galaxies to 
collide, the dissipation of angular momentum during the process 
allows the gas to be funneled to the central region of the system. 
Numerical simulations have also demonstrated how gas inflows 
associated with gas-rich or “wet” mergers can feed vigorous star 
formation (Mihos & Hernquist 1996; Chakrabarti et al. 2008). 
More recently, hydrodynamic simulations by Narayanan et al. 
(2010) uggested that SMGs can naturally form via galaxy merg¬ 
ers. From an observational point of view, this is supported by the 
clumpy or disturbed morphologies of SMGs and their complex 
kinematic signatures (e.g. Smail et al. 1998; Tacconi et al. 2008; 

Engel et al. 2010, Swinbank et al. 2011 ; Sharon et al. 2013; 

Riechers et al. 2013 ; Hung et al. 2013; Toft et al. 2014 ; 

Neri et al. 2014; Riechers et al. 2014; Chen et al. 2015). 


Engel et al. (2010) concluded that most of the SMGs with IR 
luminosities of Lir > 5 x 10 12 L© are probably major-merger 
systems [i.e. systems where the interacting galaxies have a 
mass ratio of > 1/3 (e.g. [Casey et al. 2014| )]. An alternative 
mechanism behind galaxy formation and the fuelling of their 
star formation is the accretion of gas from the intergalactic 
medium through filamentary struc tures (t he so-called cold - 
mode accretion; IKeres et al. 20051 2009; IDekel et al. 20091) . 
Cosmological hydrodynamic simulations performed by Dave 
et al. (2010) suggest that SMGs could generally obtain their 
gas reservoirs via such accretion processes (rather than through 
mergers). In these simulations, the galaxies often had complex 
morphologies and gas kinematics - signatures often inter¬ 
preted as evidence of an ongoing merger. However, as a result 
of cold-mode accretion, an extended disk-like gas structure 
undergoing rotation is also expected, and some SMGs are 
indeed found to show such signatures (e.g. Carilli et al. 20101 
|Hodge et al. 2012[ IDe Breuck et al. 20141) . Finally, we note that 
simulations suggest that during the course of their evolution, 
SMGs can exhibit properties that are reminiscent of both 
normal star-forming galaxies and vigorous starbursts (see 
[Hayward et al. 2013a[b jq For example, numerical simulations 
( [Springel & Hernquist 2005 ) have demonstrated that a disk-like 
structure can form soon after the merging of gas-r ich galaxies 
because of the rapid cooling (see also |Hopkins et al. 2009| . This 
conforms to the idea that SMGs are a heterogeneous galaxy 
population, probably caught at different stages of evolution. 

As an SMG increases its gas reservoir (through whatever 
mechanism), its central supermassive black hole (SMBH) can 
accrete increasing amounts of gas, which is driven to the nuclear 
region (e.g. lGranato et al. 20061) . Some SMGs can therefore host 
an active galactic nucleus (AGN) as revealed by deep X-ray ob¬ 
servations ( Alexander et al. 2003 i 2005; [Wang et al. 2013 ). For 
example, Chapman et al. (2005) concluded that about 20-30% of 
radio-identified SMGs display AGN activity, and radio-detected 
SMGs indeed appear to have a higher AGN fraction than the 
general SMG population ( [Wang et al. 2013) . However, most of 
the bolometric IR luminosity of SMGs is found to originate in 
star formation activity (dust-reprocessed radiation) and not gas 
accretion onto the SMBH of an AGN. An important character¬ 
istic of these accreting central black holes is that they can in¬ 
fluence the properties of their host galaxies through radiative 
and mechanical feedback. In particular, besides the exhaustion 
of the gas reservoir, AGN feedback can lead to the “quench¬ 
ing” or shut down of the star formation (e.g . | Springel et al. 2005 [ 


IDi Matteo et al. 20051 |Hopkins et al. 2006). 

Another intriguing question is the role played by SMGs 
in galaxy evolution over cosmic time. It has been found that 
SMGs are promising candidates for the progenitors of the 
most massive, passive (i.e. with little or no ongoing star for¬ 
mation) elliptical galaxies seen in the present-day universe 
(e.g . |Filly et al. \$99\ ISwinbank et al. 20061 IFu et al. 20138; 
IToft et al. 20141 |Simpson et al. 2014) . The existence of quies¬ 
cent, red massive galaxies already at z ~ 2-3 with old stel¬ 
lar populations indicates that these galaxies have experienced 
a short-lived starburst phase in their past (e.g. IRenzini 2 006; 
|Capak et al. 2008[[Coppin et al. 2010) . High-redshift (z ~ 4 - 5) 
SMGs could well represent these galaxy precursors. Besides 
their physical characteristics, the strong clustering of SMGs 


1 When a galaxy is on the main sequence, it is often said to be a “nor¬ 
mal” star-forming galaxy. In contrast, if the galaxy has a clearly en¬ 
hanced SFR with respect to its stellar mass (i.e. outlier above the main 
sequence), it is defined to be a starburst galaxy (e.g. |Magdis et al. 2011) . 
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is consistent with th is evolutionary picture (Blainetal. 2004; 
lAravena et al. 2010al) . Toft et al. (2014) found compelling ev¬ 
idence that the evolution of the giant red-and-dead ellipticals 
observed in the nearby universe, starting from z > 3 SMGs, 
goes through a transition stage manifested as compact quiescent 
galaxies at z ~ 2 - 3. 

Before the physical properties of SMGs can be studied in de¬ 
tail, the position of the source giving rise to the (sub)mm contin¬ 
uum emission must be accurately determined. The source coun¬ 
terparts at other wavelengths can only be correctly identified if 
the exact location of the FIR/(sub)mm emission is known, which 
in practice requires the analysis of FIR or (sub)mm interfero¬ 


metric observations to achieve this goal (e.g. Frayer et al. 2000 


Younger et al. 2007 

2008, 2009; Dannerbauer et al. 2008 

Aravena et al. 20101 

3 ; Smolcic et al. 2012alb; Karim et al. 2013 

Hodge et al. 2013). 

To date, however, only a few flux- 


limited SMG samples have been followed up with interfer¬ 
ometers ([Younger et al. 2007] 2009 ; ISmolcic et al. 2012bt 


[Barger et al. 2012[lKarim et al. 20131 [Hodge et al. 2013[ ). 


In this paper, we present the results of our intermediate- 
resolution (178) Plateau de Bure Interferometer (PdBI; 
IGuilloteau et al. 19921 1.3 mm continuum imaging of a sam¬ 
ple of 15 SMGs discovered by Scott et al._(2008) in the Cos¬ 
mic Evolution Survey (COSMOS; IScoville et al. 2007al) field. 
The good angular resolution of the present data allows us to 
accurately pinpoint source positions and match them with cor¬ 
rect multi wavelength counterparts. Accurate SMG positions are 
needed for their targeted spectroscopic redshift measurements, 
and knowing the proper multiwavelength counterparts allows us 
to determine the photometric redshifts of the sources - a pre¬ 
requisite for a rigorous analysis of the physical properties. After 
describing the source sample, observations, data reduction, and 
ancillary data in Sect. 2, the direct observational results and anal¬ 
ysis are presented in Sect. 3. In Sect. 4, we present our analysis of 
the redshifts of our SMGs. We then discuss our results in Sect. 5, 
and a summary is given in Sect. 6. 

In the present paper, we adopt a concordance A cold 
dark matter (ACDM) cosmology, with the Hubble constant 
Ho = 71 km s -1 Mpc -1 [i.e. the reduced Hubble constant 

h = i7 0 /(100km s _1 Mpc -1 ) = 0.71], total (dark+luminous 
baryonic) ma tter density Q m = 0.27, and dark en ergy density 
Qa = 0.73 ( [Spergel et al. 2007|ILarson et al. 2011b . In this spa¬ 
tially flat universe, 1" corresponds to a physical spatial scale 
of 8.04, 8.48, and 7.83 kpc at redshifts of z = 1,2, and 3, re¬ 
spectively. The corresponding cosmic times are 5.94, 3.34, and 
2.19 Gyr. Magnitudes in the present paper refer to the AB mag¬ 
nitude system (see lOke 1974k 


2. Observations, data, and data reduction 

2.1. Source sample 

Our new PdBI 1.3 mm observations, described in the next sub¬ 
section, were made towards the SMGs listed in Table [T] These 
SMGs were originally discovered in the 1.1 mm continuum sur¬ 
vey of a north-west subfield (0.15 deg 2 in size) of the 2 deg 2 
COSMOS field by Scott et al. (2008). The survey was carried out 
with the Aztronomical Thermal Emission Camera, or the AzTEC 
bolometer array (IWilson et al. 2008k on the 15 m James Clerk 
Maxwell Telescope (JCMT), and the target field was centred on 
a prominent large-scale structure traced by the galaxy overden¬ 
sity dScoville et al. 2007bt that includes a massive (~ 10 15 M 0 ) 
galaxy cluster at a redshift of z - 0.73 ( Guzzo et al. 2007 ). In 


Table 1 . Source list. 


Source 

JCMT/AzTEC ID 

cdb 

° 1.1mm 

[mJy] 

S/N 1 _ 1 mm 

AzTEC16 

AzTEC_J095950.29+024416.1 

3.9 ± 1.3 

4.5 

AzTEC17 

AzTEC_J095939.30+023408.0 

3.8 ± 1.4 

4.4 

AzTEC18 

AzTEC_J095943.04+023540.2 

3.s:i| 

4.3 

AzTEC19 

AzTEC_ J100028.94+023200.3 

3 - 8 - 

4.3 

AzTEC20 

AzTEC _ J100020.14+024116.0 

3 -s4 

4.3 

AzTEC21 

AzTEC-J100002.74+024645.0 

3.4:lf 

4.2 

AzTEC22 

AzTEC-J095950.69+022829.5 

3 6 +15 

J -16 

4.2 

AzTEC23 

AzTEC_ J095931.57+023601.5 

3 4 + i-4 
-i 5 

4.1 

AzTEC24 

AzTEC_J100038.72+023843.8 

o o + 1.4 

4.1 

AzTEC25 

AzTEC_J095950.41+024758.3 

3.3 + i.4 

4.1 

AzTEC26 

AzTEC_J095959.59+023818.5 

0 0 + 1.4 

j. j_i 5 

4.0 

AzTEC27 

AzTEC-J100039.12+024052.5 

0 0 + 1.4 

4.0 

AzTEC28 

AzTEC-J100004.54+023040.1 

3 3 +1 ^ 

4.0 

AzTEC29 

AzTEC-J100026.68+023753.7 

0 o + 1.4 

4.0 

AzTEC30 

AzTEC J100003.95+023253.8 

0 0 + 1.4 

J -1.6 

4.0 


Notes. The 1.1 mm flux densities listed in column (3) are deboosted 
values from Scott et al. (2008; their Table 1) and the quoted errors rep¬ 
resent the 68% confidence interval. 


total, Scott et al. (2008) reported 50 candidate SMGs with a de¬ 
tection signal-to-noise ratio S/Ni.i mm > 3.5 (see their Table 1). 

While our PdBI observations targeted the COSMOS/AzTEC 
SMGs AzTEC 16-30, the 15 brightest SMGs detected by Scott 
et al. (2008), i.e. AzTECl-15 (S/Ni.i mm > 4.6), had previ¬ 
ously been imaged (and detected) with the Submillimetre Ar¬ 
ray (SMA) at 890 pm (2" angular resolution) by Younger et al. 
(2007, 2009). Spectroscopic and/or photometric redshifts have 
been assigned to them by Younger et al. (2007, 2009), Riechers 
et al. (2010), Capak et al. (2011), Smolcic et al. (2011, 2012b), 
and M. S. Yun et al. (in prep.). Combining these data with the 
present observations towards AzTEC 16-30 provides us with a 
flux-limited sample of 30 SMGs (S/Ni.i mm > 4), which have 
all been followed up with intermediate-resolution interferomet¬ 
ric observations. This allows us to carry out a statistically mean¬ 
ingful study of their redshift distribution. 

2.2. Intermediate-resolution 1.3 mm imaging 

The PdBI 1.3 mm (230.5 GHz) continuum observations towards 
AzTEC 16-30 (project W0AE) were carried out between January 
and November 2013. The array of six antennas was mostly in 
its C configuration, i.e. the second most-compact configuration 
with 15 baselines ranging from 24 to 176 m (which corresponds 
to 18.5-135.4 kT). On 16 April, however, antenna station E04 
was not available, resulting in only ten baselines. On 12 October, 
when part of the observations towards AzTEC 16-22 were per¬ 
formed, the array was in its most compact D configuration (base¬ 
lines in the range 24-97 m or 18.5-74.6 k A). The lower-sideband 
(LSB) system temperature was typically ~ 200 K, except on 
16 April and 3 May when it was 300-350 K. The atmospheric 
precipitable water vapour (PWV) was typically in the range 2- 
4 mm, except on 16 April when it was 6 mm. The best weather 
conditions were on 30 October, when the PWV value was only 
about 1 mm. The phase centres used were the AzTEC 1.1 mm 
peak positions of the sources detected by Scott et al. (2008), and 
the on-source observing time per source was ~ 1.5 hr. 

The backend used was the WideX correlator, which is com¬ 
posed of four units [two for both orthogonal linear (horizontal 
and vertical) polarisation modes], each being 2 GHz wide and 
split into 1 024 channels (corresponding to a channel width of 
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about 1.95 MHz). The total effective bandwidth is about 1.8 GHz 
for each unit or about 3.6 GHz for both polarisations. The cor¬ 
relator bandpass calibration was achieved by observing 3C84 
(NGC 1275), 3C279, and B0923+392. Gain phases and ampli¬ 
tudes were calibrated by observing B0906+015 and 1005+066. 
The standard source 3C84 was used the most often as a flux cal¬ 
ibrator, with the adopted 230 GHz flux density of S 230 GHz = 
9.85 - 12.36 Jy depending on the observing day. The other flux 
calibrators used were 3C279 (10.68 Jy), 0851+202 (3.04 Jy), 
and 0923+392 (2.5-3.16 Jy). The absolute flux-calibration un¬ 
certainty was estimated to be about 20%, which is based on the 
observed scatter in the calibrators’ flux densities. The primary 
beam (PB) of the PdBI at the observing frequency is 2173 (Half- 
Power Beam Width or HPBW). At z = 2, this corresponds to 
about 180.5 physical kpc. 

Calibration and imaging were performed using the CLIC 
(Continuum and Line Interferometer Calibration) and MAP¬ 
PING programs of the GILDAS software packag^] respectively. 
When creating the maps, natural weighting was applied to the 
calibrated visibilities (i.e. weighting according to the number 
of measurements within a given region of the wv-plane). The 
CLEAN algorithm was used for deconvolution, and applied in 
regions centred on the strongest emission features. The typi¬ 
cal resulting synthesised beam size (Full Width at Half Maxi¬ 
mum or FWHM) is 178, and the restored continuum maps (dual 
polarisation) have l<x root mean square (rms) noise values of 
~ 0.2 mJy beam -1 . Hence, the statistical positional error, A# stat ^ 


a/# majflmin/(2S/N) (where flmaj and 0 m i n are the major and minor 
axes of the beam; Rei d et al. 19 88; Younger et al. 2007), can be 
estimated to be about 0718/S' i. 3 mm[inJy beam -1 ], which is about 
0723 for a typical 4 cr source. We note that merging the C and 
D configuration visibilities together for AzTEC 16-22 improved 
the uv coverage and produced maps of higher spatial dynamic- 
range than those of AzTEC23-30. The phase centres, both the 
synthesised beam sizes and position angles, and the rms noises 
of the final cleaned, continuum maps are listed in Table [2] 


2.3. The COSMOS field - ancillary data 

Since our target sources lie within the COSMOS field, they have 
been observed with several ground- and space-based telescopes 
at wavelengths spanning from the X-rays to the radio regime. 

Observations at X-ray wavelengths were performed with the 
XMM-Newton ( Hasing er et al. 2067]), and Chandra satellites (C- 
COSMOS; [Elvis et al. 2009t IPuccetti et al. 20091) . The Galaxy 
Evolution Explorer , or GALEX , imaged the COSMOS field in 
the near-UV (NUV) and far-UV (FUV) ( [Zamojski et al. 2007[ ). 
Hubble Space Telescope (/AST) Advanced Camera for Sur¬ 
veys (ACS) Wide Field Channel (WFC) observations of the 
COSMOS field were done in the I band (the F814W filter) 
dScoville et al. 2007al IKoekemoer et al. 20071) . Observations at 
optical/near-IR (NIR) wavelengths have been carried out with 
the 8.2 m Subaru telescope, the 3.6 m Canada France Hawaii 
Telescope (CFHT), the 3.8 m United Kingdom Infrared Tele¬ 
scope (UKIRT), the 2.2 m University of Hawaii telescope called 
UH88 (or UH2.2), and the 4 m telescopes of the Kitt-Peak Na¬ 
tional Observatory (KPNO), the Cerro Tololo Inter-American 
Observatory (CTIO), and the National Optical Astronomy Ob¬ 
servatory (NOAO) [see Capak et al. (2007), Taniguchi et al. 
(2007), and McCracken et al. (2010) for details]. 


2 Grenoble Image and Line Data Analysis Software is pro¬ 
vided and actively developed by IRAM, and is available at 
http://www.iram.fr/IRAMFR/GILDAS 


Table 2. Observational parameters. 


Source 

<* *2000.0 

[h:m:s] 

£2000.0 

$syn 

["x"] 

P.A. 

n 

<* rms 

[mJy beam -1 ] 

AzTEC16 

09 59 50.29 

+02 44 16.1 

1.87 x 1.18 

12.82 

0.247 

AzTEC17 

09 59 39.30 

+02 34 08.0 

1.85 x 1.15 

13.49 

0.239 

AzTEC18 

09 59 43.04 

+02 35 40.2 

1.85 x 1.14 

13.75 

0.256 

AzTEC19 

10 00 28.94 

+02 32 00.3 

1.85 x 1.14 

13.67 

0.302 

AzTEC20 

10 00 20.14 

+02 41 16.0 

1.85 x 1.14 

14.42 

0.252 

AzTEC21 

10 00 02.74 

+02 46 45.0 

1.86 x 1.14 

14.82 

0.256 

AzTEC22 

09 59 50.69 

+02 28 29.5 

1.80 x 1.21 

0.00 

0.227 

AzTEC23 

09 59 31.57 

+02 36 01.5 

1.76 x 1.03 

4.19 

0.205 

AzTEC24 

10 00 38.72 

+02 38 43.8 

1.76 x 1.03 

4.53 

0.188 

AzTEC25 

09 59 50.41 

+02 47 58.3 

1.75 x 1.03 

5.00 

0.191 

AzTEC26 

09 59 59.59 

+02 38 18.5 

1.76 x 1.03 

5.29 

0.178 

AzTEC27 

10 00 39.12 

+02 40 52.5 

1.76 x 1.02 

5.19 

0.215 

AzTEC28 

10 00 04.54 

+02 30 40.1 

1.76 x 1.02 

5.28 

0.225 

AzTEC29 

10 00 26.68 

+02 37 53.7 

1.75 x 1.02 

5.57 

0.212 

AzTEC30 

10 00 03.95 

+02 32 53.8 

1.77 x 1.01 

6.09 

0.205 


Notes. The equatorial coordinates refer to the JCMT/AzTEC 1.1 mm 
peak positions ( Scott et al. 2008 ). and they represent the phase centres 
of the PdBI observations presented here. 


New NIR imaging of the COSMOS field in the Y (1.02 pm), 
J (1.25 pm), H (1.65 pm), and K s (2.15 pm) bands is be¬ 
ing collected by the UltraVISTA survey ([McCracken et al. 20121 
lllbert et al. 20130 Mid-infrared (MIR) observations were ob¬ 
tained with the Infrared Array Camera (IRAC; 3.6-8.0 pm) 
and the Multiband Imaging Photometer for Spitzer (MIPS; 24- 
160 pm) on board the Spitzer Space Telescope as part of the 
COSMOS Spitzer survey (S-COSMOS; Sanders et al. 2007). 
Far-infrared (70, 160, and 250 pm) to submm (350 and 
500 pm) Herschel continuum observations were performed 
as part of the Photodetector Array Camera and Spectrom¬ 
eter (PACS) Evolutionary Probe (PEP; [Lutz et al. 20111 and 
the Herschel Multi-tiered Extragalactic Survey (HerMES0 
IQliver et al. 2012b programmes. Radio-continuum imaging was 
done at 20 cm (1.4 GHz) with the Very Large Array or VLA 
(ISchinnerer et al. 20071 2010), and at 10 cm (3 GHz) with the 
Karl G. Jansky Very Large Array (VLA-COSMOS 3 GHz Large 
Project; PI: V. Smolcic). In addition to the imaging observations, 
a large spectroscopic redshift survey of galaxies in the COS¬ 
MOS field has been carried out with the Very Large Telescope 
(VLT), a survey known as zCOSMOS ( [Lilly et al. 20071 2009), 
and with the Keck telescope (M. Salvato et al., in prep.). Photo¬ 
metric redshifts towards sources in the COSMOS field have been 
computed using 30 wavelength bands spanning from UV to MIR 
(lllbert et al. 20091 2013). 

Most of the extensive multiwavelength datasets are publicly 
available from the NASA/IPAC Infrared Science Archiv^] 

3. Source catalogue, multiwavelength counterparts, 
and multiplicity 

3.1. Source extraction and multiwavelength counterparts 

The PdBI 1.3 mm images towards AzTEC 16-30 are shown in 
Fig. [I] We note that the PdBI PB at the observed frequency, 
2173, closely resembles the size of the JCMT/AzTEC beam 
of about 18" (FWHM), and that our observation wavelength 
(1.3 mm) is very close to that of the original discovery observa¬ 
tions (JCMT/AzTEC) by Scott et al. (2008; 1.1 mm), facilitating 
comparison between these two studies. 

3 The data products are produced by TERAPIX; see 
http://terapix.iap.fr 

* http://hermes.sussex.ac.uk 

5 http://irsa.ipac.caltech.edu/Missions/cosmos.html 
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To systematically search for sources in the final, cleaned 
maps, we followed Hodge et al. (2013) and Karim et al. (2013). 
Briefly, we used an Interactive Data Language (IDL)-based 
source-extraction routine that first finds pixel values above 2.5<x, 
where cr is the rms noise determined using non-overlapping rect¬ 
angular apertures across the map. The size of each aperture was 
taken to be large compared to the synthesised beam so that 
each of them will yield a representative sampling of indepen¬ 
dent beams. Since some apertures contained physical signal (i.e. 
sources), the value of cr was taken to be the median of all differ¬ 
ent rms values. The value of cr derived this way is consistent with 
the GILDAS-derived map rms noise given in Col. (6) in Table [2j 
the first value was found to be 6% higher on average than the sec¬ 
ond value (the median ratio between the two rms values is 1.05). 
The routine then models the source emission within a 3" x 3" 
region using a Gaussian and the MCMC (Metropolis-Hastings 
Markov chain Monte Carlo) algorithm. Extended sources are fit 
for six parameters (peak surface brightness, peak position, mi¬ 
nor axis, major-to-minor axis ratio, and position angle) while 
for point sources the size is fixed to that of the synthesised 
beam (leaving only three free parameters). To generate a robust 
catalogue of PdBI-detected sources we adopt the approach al¬ 
ready used by Smolcic et al. (2012b) and consider sources with 
S/N > 4.5 in the PdBI 1.3 mm maps as robust detections, while 
sources with 4 < S/N < 4.5 are considered real only if they 
have a multi wavelength counterpart within a search radius of 1'.'5 
(within the COSMOS UltraVISTA NIR, Spitzer IRAC/MIPS, 
VLA 10 or 20 cm radio catalogues). When multiple PdBI source 
candidates are detected, we label them AzTEC17a, AzTEC17b, 
etc. 


In total we find 22 sources, 11 of which are associated with 
multiwavelength counterparts (three additional sources have a 
nearby ACS I -band source but no “counterparts” at other wave¬ 
lengths). We report their positions (J2000.0 equatorial coordi¬ 
nates and offset from the phase centre) and primary-beam cor¬ 
rected flux densities in Table [3] For the unresolved sources, we 
report the peak flux density, while for the (marginally) resolved 
sources, AzTEC21a, 27, and 28, we give the total flux density 
derived from the best-fit six-parameter model. The latter val¬ 
ues were also independently derived by summing over all pixels 
within the 2.5cr contour of 1.3 mm emission, and were found to 
be consistent with the model values. The quoted flux density un¬ 
certainties are based on the rms noise values and the 20% abso¬ 
lute calibration error (added in quadrature). We note that inspec¬ 
tion of the residual maps of the model Gaussian fits showed that 
AzTEC21a could be well-fitted with a deconvolved FWHM size 
of 0majX#min = (276+172)x(073±075). However, AzTEC27 and 
28 are not as well modelled by a single Gaussian source model. 
For AzTEC27, only the major axis of the elliptical Gaussian 
could be determined (0 ma j = 376), while the size of AzTEC28 
was determined to be 0 ma j x 0 min = (175 ± 078) x (076 ± 077). 
The large uncertainties in the sizes reflect that these sources are 
only marginally resolved and/or are not well-represented by a 
single Gaussian (see e.g. [Condon 19971) . In all the three cases 
(AzTEC21a, 27, and 28), however, the peak flux density was 
found to be lower than the total flux density, supporting their 
marginally extended nature. The multi wavelength counterpart 
IDs of our PdBI sources are reported in Col. (7) in Table [3] and 
the last column lists the projected angular offset from the PdBI 
source. Multi wavelength zoomed-in images towards our sources 
are provided in Fig. |A.1| Notes on individual sources are given 
in Appendix C. 


3.2. Source catalogue validation 

To test the robustness of our PdBI source identifications we 
quantified the number of expected spurious sources by search¬ 
ing for detections in the negative part of the PdBI maps in the 
same way as described above. We found no spurious sources as¬ 
sociated with multiwavelength counterparts within a search ra¬ 
dius of 175. This is consistent with the random match probability 
within this radius, based on the optical/IR/radio catalogue source 
densities, of only ~ 0.2% (tSmolcic et al. 2012bl) . It also implies 
that all our PdBI detections with multi wavelength counterparts 
are likely to be real. We also find that the number of spurious 
sources increases with increasing distance from the phase cen¬ 
tre, as shown in Fig. [2] Out to a distance of 6" from the phase 
centre only one spurious source is expected. This suggests that 
potentially one of the three sources (AzTEC20, 26b, 28; Table[3]) 
detected with S/N > 4.5, within 6" from the phase centre and 
with no multi wavelength counterparts may be spurious. At a 
distance of 9 - 13" from the phase centre a total of five spu¬ 
rious sources with S/N > 4.5 is expected. This suggests that 
five sources (AzTEC22, 24abc, 30; see Table [3]) we detect at 
S/N > 4.5 and within 973 - 1278 from the phase centre and 
with no multi wavelength counterparts may be spurious. 

To test the validity of our sources further, we compare for 
each PdBI source the PdBI 1.3 mm flux to the AzTEC/JCMT 
1.1 mm flux in Fig. [3] For this, the deboosted AzTEC 1.1 mm 
flux densities given in Table [T] were scaled down using the com¬ 
mon assumption that the dust emissivity index is J3 = 1.5 (e.g. 
IDunne & Eales 20011 [Barger et al. 2012| and references therein). 
In case multiple PdBI sources are extracted from the target 
field, the sum of their flux densities is plotted (AzTEC25 was 
not detected in the PdBI map). As can be seen, the two val¬ 
ues are generally in reasonable agreement with each other. We 
note, however, that a one-to-one correspondence is not expected 
given observational limitations (such as flux deboosting methods 
and possible blending of multiple sources in the low-resolution 
single-dish data; see e.g. Fig. 2 in ISmolcic et al. 2012bl for a 
comparison of SMA and JCMT/AzTEC 890 pm flux densi¬ 
ties for AzTEC 1-15). The four sources showing higher PdBI 
flux densities than expected from single-dish measurements are 
AzTEC19, 21, 24, and 29. Given that AzTEC24 has no multi¬ 
wavelength counterpart associated, and that its PdBI flux is in¬ 
consistent with that from the JCMT/AzTEC survey may sug¬ 
gest that the three identified PdBI sources within the AzTEC24 
field are spurious. On the other hand, emission features at 
the border of and/or just outside the PB FWHM (as in the 
cases of AzTEC 19, 24, and 29) may not be contributing to the 
JCMT/AzTEC source flux densities detected at 18" resolution, 
and deboosted JCMT/AzTEC flux densities may also be invalid 
in this comparison for sources with widely separated compo¬ 
nents. In summary, we conclude that 4 out of the 22 PdBI iden¬ 
tified sources may be spurious. 

3.3. Multiplicity of single-dish detected, JCMT/AzTEC SMGs 

Submillimetre sources identified in single-dish studies can be 
composed of multiple components, and this multiplicity can 
be revealed by higher-resolution (interferometric) observations. 
These components can typically be associated with individual 
galaxies that might be physically related (and potentially inter¬ 
acting), or might be just chance alignments of galaxies lying at 
different redshifts. As the multiplicity fraction of submm sources 
depends on the initial resolution of the single-dish observations 
and on the depth of the interferometric follow-up, it is not sen- 
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big. 1 . The Tctbl r.5 mm^ehtinuum images towards AzTEC16-30, annotated with the source designations (AzTEC17a, 17b, etc.). The images 
are shown with linear scaling, and the overlaid solid contours start from 3cr and increase in steps of lcr. The white dashed contours show the 
corresponding negative features (starting from -3cr). Each image is centred on the phase centre position, i.e. the original AzTEC 1.1 mm centroid 
from Scott et al. (2008) marked by the plus sign. The dashed circle shows the primary beam FWHM of 2173, while the filled ellipse in the bottom 
left corner shows the synthesised beam size. The images are not corrected for the primary beam attenuation. 
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sible to provide a simple definition for multiplicity. The present 
study benefits from the fact that the PB of our PdBI follow-up 
observations is well matched to the beam FWHM of the JCMT 
observations used for identifying the initial source sample (21'.'3 
versus 18"; Sect. 3.1). This facilitates the analysis of source 
blending in the single-dish maps. Based on our PdBI source 
detections six out of the 15 AzTEC 16-30 SMGs (or 5/15 if 
AzTEC24 is treated as spurious) separate into multiple compo¬ 
nents. This corresponds to 40 ± 16% (33 ± 15%), where the er¬ 
rors are Poisson errors. Within the sample of the brightest COS¬ 
MOS JCMT/AzTEC SMGs (AzTECl-15) only 2/15 (13 ± 9%) 
separate into multipl e components ( [Younger et al. 2007] 2009; 
ISmolcic et al. 2012bl) . Combining these two samples then yields 
that 7-8 SMGs in our flux- and S/N-limited single-dish de¬ 
tected, COSMOS JCMT/AzTEC sample of 30 sources separate 
into multiple sources when observed at < 2" angular resolution. 
This corresponds to ~ 25 ± 9%, and will be discussed further in 
Sect. 5.1. 



Fig. 2. Cumulative number distribution of spurious sources identified 
in the negative parts of the PdBI maps (S/N < -4.5) as a function 
of angular distance ["] from the phase centre. The vertical dashed line 
marks the half-power radius of the PB. 


4. Redshift distribution of 1.1-mm selected SMGs in 
the COSMOS field 

4.1. Redshifts for AzTEC1-15 

Among the SMA-detected SMGs AzTEC 1-15, there are eight 
spectroscopic redshifts reported in the literature (AzTEC 1, 2, 3, 

5, 6, 8, 9, and 11; see lSmolcic et al. 2012bt their Tables 1 and 4 
and references therein). However, as described in Appendix B, 
the spectroscopically determined redshifts for AzTEC5, 6, and 
9 are uncertain because of the poor quality of the spectrum or 
contamination by foreground galaxies. Besides these cases, we 
discuss the updated redshifts among AzTEC 1-15 in Appendix B 
(the redshifts are listed in Table]?]). 

4.2. Redshifts for AzTECI6-30 

The optical/IR photometric redshifts of AzTEC 16-30 were com¬ 
puted (when possible) by fitting optimised spectral templates to 
their spectral energy distributions (SEDs) using the HYPERZ code 



16 18 20 22 24 26 28 30 

AzTEC SMG ID 


Fig. 3. A plot showing the PdBI/JCMT 1.3 mm flux density ratio for 
the SMGs AzTEC 16-30. The single-dish JCMT/AzTEC flux densities 
were scaled from the deboosted 1.1 mm flux densities (Table [TJ by as¬ 
suming that the dust emissivity index is p = 1.5. The dashed line repre¬ 
sents the line of equality. 


dBolzonella et al. 2000| ^] The redshift was treated as a free pa¬ 
rameter and determined using a x 2 minimisation method, i.e. 
the most likely redshift was determined statistically (see also 
ISmolcic et al. 2012al b). We used the Calzetti et al. (2000) ex¬ 
tinction law, and the V-band extinction, Ay, was varied from 0 
to 5 mag. The allowed redshift range was z e [0, 7]. The spec¬ 
tral templates used were generated with the GRASIL radiative- 
transfer code dSilva et al. 19981 [Tglesias-Paramo et al. 2007 } and 
optimised for SMGs using the method described by Michalowski 
et al. (2010). When comparing results for their tested spectral 
templates, Smolcic et al. (2012a) found that the best agreement 
(i.e. the tightest x 2 distribution) between the spectroscopic and 
photometric redshifts (sample of eight COSMOS SMGs) is ob¬ 
tained when employing the templates derived by Michalowski 
et al. (2010). Smolcic et al. (2012b; see their Fig. 7) repeated 
the analysis using a larger source sample, and their similar re¬ 
sults lend further support to the reliability of the Michalowski 
et al. (2010) spectral-template library. We therefore decided to 
perform our photo-z analysis using this library of templates. The 
optical/IR SEDs for the identified SMG counterparts are shown 
in Fig. [4] The template-fitting method of finding the best photo-z 
solution is based on the minimisation of the reduced chi-square 
(Yred) va l ue > which is the chi-square divided by the number of 
degrees of freedom (dof) [see Eq. (1) in lBolzonella et al. 2000M . 
The HYPERZ program yields the probability associated with the 
minimum x 2 ed for each redshift step, P(z) = exp(-^ ed /2). The 
absolute (total) chi-square (XtoO distribution for each source as 
a function of redshift is presented in a panel next to the corre¬ 
sponding SED plot in Fig. [4] We computed the formal lower and 
upper 99% confidence limits for the best-fit photo-z value. For¬ 
mally, these were calculated from the x 2 probability distribution 
P(X 2 < \ 2 \v) = 0.99 ( Avni 1976 ). where A x 2 is the increment 
in x 2 required to cover the parameter space region with a z con¬ 
fidence of 99%, and v is the number of dof. The confidence in¬ 
terval (Cl) equals the set of all z values that satisfy the condition 
X 2 (Z) -Xlin ^ V- 

For those sources with no optical/IR counterparts, the above 
method could not be used to derive their redshift. The sources 


6 http://webast.ast.obs-mip.fr/hyperz/ 
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Table 3. Characteristics and optical-to-MIR and VLA 20 cm counterparts of the SMGs identified in the PdBI 1.3 mm maps. 


Source 

2000.0 

[hums] 

£2000.0 

[V:"] 

S 1.3 mm 

[mJy] 

S/Nomm 

Offset 

["] 

Candidate 
counterpart ID a 

r 

["] 

AzTEC 16 

09 59 50.069 

+02 44 24.50 

2.07 + 0.62 

5.0 

9.0 

501 (ACS-/) 

0.35 

AzTEC 17a 

09 59 39.194 

+02 34 03.83 

1.58 + 0.43 

6.2 

4.5 

613229 (1 GALEX) 

1.29 







1496 (ACS-/) 

1.36 







1475165 (COSMOS+UltraVISTA) 

1.38 







271694 (New UltraVISTA) 

1.47 







250117 (IRAC) 

0.30 







J095939.19+023403.6 (VLA Deep) 

0.22 

AzTEC 17b 

09 59 38.904 

+02 34 04.69 

1.53 + 0.46 

4.5 

6.8 

1475223 (COSMOS+UltraVISTA) 

1.25 

AzTEC 18 

09 59 42.607 

+02 35 36.96 

1.78 + 0.54 

4.5 

7.3 

1044 (ACS-/) 

0.76 







1471053 (COSMOS+UltraVISTA) 

0.82 

AzTEC 19a 

10 00 28.735 

+02 32 03.84 

3.98 + 0.91 

10.3 

4.4 

1593 (ACS-/) 

0.40 







1455882 (COSMOS+UltraVISTA) 

0.20 







262214 (New UltraVISTA) 

0.20 







242501 (IRAC) 

0.63 







2158 (MIPS 24 pm) 

0.15 







J100028.70+023203.7 (VLA Deep) 

0.53 

AzTEC 19b 

10 00 29.256 

+02 32 09.82 

5.21 + 1.30 

9.7 

10.6 

597821 0 GALEX) 

0.66 







1486 (ACS-/) 

0.68 







1455681 (COSMOS+Ultra VISTA) 

0.63 







262766 (New UltraVISTA) 

0.74 

AzTEC20 

10 00 20.251 

+02 41 21.66 

1.85 + 0.49 

6.0 

5.9 



AzTEC21a 

10 00 02.558 

+02 46 41.74 

3.37 + 1.03 

5.8 

4.3 

711447 {GALEX) 

0.35 







1688587 (COSMOS+UltraVISTA) 

0.15 







328878 (New UltraVISTA) 

1.00 







297396 (IRAC) 

1.11 







7262 (MIPS 24 pm) 

1.46 

AzTEC21b 

10 00 02.710 

+02 46 44.51 

1.34 + 0.38 

4.2 

0.7 

711786 {GALEX) 

0.53 







1688585 (COSMOS+UltraVISTA) 

0.45 

AzTEC21c 

10 00 02.856 

+02 46 40.80 

1.27 + 0.40 

4.5 

4.5 

712026 {GALEX) 

1.49 







297223 (IRAC) 

1.46 

AzTEC22 

09 59 50.681 

+02 28 19.06 

1.82 + 0.59 

5.1 

10.5 



AzTEC23 

09 59 31.399 

+02 36 04.61 

0.99 + 0.29 

4.9 

4.0 

1494354 (COSMOS+UltraVISTA) 

1.24 

AzTEC24a 

10 00 38.969 

+02 38 33.90 

1.79 + 0.53 

4.9 

10.6 



AzTEC24b 

10 00 39.410 

+02 38 46.97 

1.72 + 0.53 

5.0 

10.8 



AzTEC24c 

10 00 39.194 

+02 38 54.46 

2.85 + 0.78 

5.1 

12.8 



AzTEC25 z? 








AzTEC26a 

09 59 59.386 

+02 38 15.36 

0.98 + 0.28 

5.4 

4.4 

647670 {GALEX) 

1.05 







930 (ACS-/) 

1.07 







1709726 (COSMOS+UltraVISTA) 

0.94 







291786 (New UltraVISTA) 

0.83 

AzTEC26b 

09 59 59.657 

+02 38 21.08 

0.90 + 0.26 

4.8 

2.8 



AzTEC27 

10 00 39.211 

+02 40 52.18 

3.36 + 0.97 c 

6.0 

1.4 

666 (ACS-/) 

1.15 

AzTEC28 

10 00 04.680 

+02 30 37.30 

2.38 + 0.77 

5.5 

3.5 



AzTEC29a 

10 00 26.351 

+02 37 44.15 

2.45 + 0.67 

4.7 

10.8 

736 (ACS-/) 

1.41 

AzTEC29b 

10 00 26.561 

+02 38 05.14 

9.01+2.39 

7.3 

11.6 

1685295 (COSMOS+UltraVISTA) 

0.76 

AzTEC30 

10 00 03.552 

+02 33 00.94 

1.53 + 0.45 

4.6 

9.3 




Notes. The coordinates given in columns (2) and (3) refer to the peak position determined by the three-parameter point-source model fit for 
all sources with peak surface brightness of S/N > 4. The flux densities given in column (4) are primary-beam corrected, and the quoted errors 
include the 20% calibration uncertainty. For the (marginally) resolved sources AzTEC21a, AzTEC27, and AzTEC28, the total flux density 
was derived from the best-fit six-parameter source model. The S/N ratio given in column (5) refers to the extracted value in the non-primary- 
beam-corrected map. In column (6), we give the PdBI source offset from the phase centre, i.e. the AzTEC 1.1 mm centroid. The last column 
gives the projected angular separation between the 1.3 mm peak position and the counterpart position. (a) The references for the different source 
catalogues are as follows: ACS I -band dLeauthaud et al. 20071 ): GALEX (M. A. Zamojski et al., in prep.); COSMOS photometry catalogue 
UCapak et al. 2007]); UltraVIS TA ( |Capak et al . 2007 McCracken et al. 2012llllbert et al. 20131) : Spitzer IRAC/MIPS (S-COSMOS team); VLA 
Deep dSchinnerer et al. 20101 ). (Z?) The 1.3 mm features in AzTEC25 did not fulfil our detection criteria. (c) AzTEC27 is probably subject to 
gravitational lensing (see Appendices C and D), and our lens modelling suggests a magnification factor of p - 2.04 + 0.16. In this case, 
AzTEC27’s intrinsic flux density at observing-frame 1.3 mm would be 1.65 + 0.49 mJy. 


that are not spurious are likely to lie at high redshift or/and be 
heavily obscured by dust. Since the radio and submm contin¬ 
uum have very different ^-corrections, the radio/submm flux- 
density ratio strongly depends on the source redshift. As pro¬ 
posed by Carilli & Yun (1999, 2000), the spectral index between 
860 pm (350 GHz) and 20 cm (1.4 GHz), c^ 5 4 °, can be used to 
estimate the redshift. The 860 pm flux densities were estimated 
from the 1.3 mm values by assuming that the dust emissivity 
index is (3 = 1.5. On the basis of this, we used the mean radio- 
to-submm spectral index, (or| 5 4 0 ), predicted from 17 low-redshift 
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star-forming galaxies by Carilli & Yun (2000), to set constraints 
on the source redshift. For 1.4 GHz non-detections, we used the 
3<x upper limit to 5 i.4ghz (typically < 0.05 mJy) to derive a 
lower limit to a 35 ®, hence a lower limit to z. The uncertainty in 
the radio/submm redshift was determined from those associated 
with the flux densities. For AzTEC27, which is subject to gravi¬ 
tational lensing, the differential lensing effects were assumed to 
be negligible (i.e. the boost in flux density was assumed to be in¬ 
dependent of wavelength), hence the radio/submm flux density 
ratio was assumed to be independent of the magnification factor. 
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We note that a 1.4 GHz non-detection could simply be related to 
a spurious PdBI 1.3 mm source. Therefore, some of the derived 
lower limits to z should be taken with caution. 

The redshifts of AzTEC 16-30 are discussed in more detail in 
Appendix C. In summary, for one source (AzTEC 17a) we have 
a spectroscopic redshift, for nine sources we have derived pho¬ 
tometric redshifts, and for 12 sources we have submm-radio flux 
ratio based redshift estimates. 

4.3. Redshift distribution of AzTECI-30 

All the derived and adopted redshifts for AzTEC 1-30 are listed 
in Table [4] n total, for six (five among AzTEC 1-15, one among 
AzTEC 16-30) sources in the sample we have a spectroscopic 
redshift, for 17 a photo-z (eight among AzTEC 1-15, nine among 
AzTEC 16-30), while for 15 (three for AzTECl-15, 12 for 
AzTEC 16-30) we have a redshift estimated from the submm- 
radio flux-density ratio. By combining the up-dated redshifts 
of AzTEC 1-15 with the present results, we derived the red¬ 
shift distribution for the combined sample of AzTEC 1-30. The 
constructed redshift distributions are shown in the top panel of 
Fig. [5] The bottom panel of Fig. [5] shows the probability density 
distribution [P oc exp(-^ 2 /2); a kernel density estimate], of our 
total SMG sample constructed using the same redshift data as 
in the histograms in the top panel. The probability distribution 
functions (PDFs) were summed for i) sources with z spec values, 
where the PDF was assumed to be a delta function [S = c>(z sp ec)]> 

ii) sources with z p hot values (PDFs derived using HYPERZ), and 

iii) sources that had only lower limits to z, in which case the PDF 
was assumed to be a flat function from the lower z limit up to 
z = 6.5. Before calculating the overall PDF, the individual PDFs 
were normalised so that their integral becomes unity. From this 
distribution we derived a median redshift of z = 3.20 and 68% 
Cl of z = 1.39-5.26. 

We further calculated the statistical parameters indepen¬ 
dently using the R program package called Nondetects And Data 
Analysis for environmental data (NADAi fHelsel 20051) . which is 
an implementation of the statistical methods provided by the As¬ 
tronomy Survival Analysis (ASURV; [Feigelson & Nelson 1985| 
package. This method robustly takes lower redshift limits into 
account (e.g. Yun et al. 2012 ). It was assumed that the censored 
data follow the same distribution as the actual values, and we 
used the Kaplan-Meier (K-M) method to construct a model of 
the input data. After applying the K-M survival estimator, we 
found that the mean, median, standard deviation (std), and 95% 
Cl of the redshifts of AzTECl-15 are <z> = 3.16 ± 0.37, z = 
3.05 ± 0.44, std = 1.48, and Cl = 2.44 - 3.89, respectively. For 
AzTEC16-30, the corresponding values are (z) = 3.02 ± 0.20, 
z = 3.20 ± 0.25, std = 0.92, and Cl = 2.63 - 3.40. The 
combined sample (open grey histrogram in Fig. [5} has the val¬ 
ues (z) = 3.19 ± 0.22, z = 3.17 ± 0.27, std = 1.35, and 
Cl = 2.76 - 3.62. We note that the median redshift of 3.20 de¬ 
rived from the redshift PDF shown in the bottom panel of Fig. [5] 
is in excellent agreement with the value we derived using the 
survival analysis. The median redshift, z = 3.17 ± 0.27, corre¬ 
sponds to a luminosity distance of Jl = 27.6^7 Gpc. Finally, 
we performed a two-sided Kolmogorov-Smirnov (K-S) test be¬ 
tween the z distributions of AzTEC 1-15 and AzTEC 16-30. By 
excluding the lower limits from the sample^] the maximum dif¬ 
ference between the cumulative distributions, i.e. the D statistic, 
was found to be D = 0.2615, while the associated probability 

7 In the survival analysis it was assumed that the lower limits and exact 
values have a common distribution. 


that the two distributions are drawn from the same parent distri¬ 
bution is p - 73.4%. Moreover, the Welch’s two-sample Mest 
under the null hypothesis that the two means are equal yields a 
/7-value of about 0.338 (when excluding the lower limits), mean¬ 
ing that there is no evidence for a difference in the mean redshifts 
of AzTEC 1-15 and AzTEC 16-30. However, as shown in Fig. [5] 
the highest redshift SMGs (z > 4.3) in our sample are found 
among AzTEC 1-15. The redshift distributions derived in other 
SMG surveys, and how they compare to the present results, will 
be discussed in Sect. 5.3. 


5. Discussion 

5.1. PdBi 1.3 mm imaging results and source multiplicity 

Even though our SMGs were detected at 4 - 4.5<x significance 
in the JCMT/AzTEC 1.1 mm survey dScott et al. 2008k not all 
of them are (clearly) detected in the present higher-resolution 
PdBi 1.3 mm imaging. It is possible that some of the weak/non- 
detected AzTEC sources are actually composed of multiple ob¬ 
jects that are too faint to be detected at the current detec¬ 
tion limit [see Large APEX BOlometer CAmera (LAB OCA) 
compared to Atacama Large Millimetre/submillimetre Array 
(ALMA ); [Hodge et al. 2013[ and SCUBA compared to SMA; 
IChen et al. 2014M . As can be seen in Fig. [3j only one source 
(AzTEC23) appears to have < 1, so we are not 

generally missing extended emission in our SMGs. Another rea¬ 
son for some of our PdBi non-detections could be that some of 
the original JCMT/AzTEC detections are spurious. Out of the 50 
JCMT/AzTEC SMGs rerported by Scott et al. (2008), 48 (96%) 
lie within the region mapped with AzTEC on the 10 m Atacama 
Submillimetre Telescope Experiment (ASTE) by Aretxaga et al. 
(2011, 2012). However, only 16 JCMT/AzTEC-detected sources 
are common to the ASTE/AzTEC 1.1 mm catalogue (Table 1 in 
[Aretxaga et al. 20lT| . For example, the ASTE/AzTEC 1.1 mm 
image extracted from the position of AzTEC30 shows emis¬ 
sion only in the eastern part of the target region. Although the 
difference is, at least partly, caused by the difference in angu¬ 
lar resolutions of ASTE (34"; [Aretxaga et al. 2011| ) and JCMT 


(18"), it is possible that the JCMT/AzTEC source is spurious. 
The expected false-detection rate in the JCMT/AzTEC survey at 
S/Ni.imm > 4.0 is < 2 sources (IScott et al. 20081 Fig. 7 therein). 
Because all the 15 brightest sources AzTEC 1-15 are interfero- 
metrically confirmed ( [Younger et al. 2007 [ 2009), we could ex¬ 
pect to find about two spurious sources among AzTEC 16-30. 

When observed with single-dish telescopes, multicomponent 
source systems can be blended and give an impression of a sin¬ 
gle source. This can be the case even if the sources are not 
physically related to each other, i.e. they can lie at significantly 
different redshifts (e.g. Cowley et al. 201 5| . Among our target 
fields, AzTEC 17, 19, 21, 24, 26, and 29 appear to show two to 
three source components. Of the 15 observed AzTEC single-dish 
detected SMGs, this would mean that 40 ± 16% are multiple 
systems (or, as explained in Sect. 3.3, 33 ± 15% if AzTEC24 
is not included)]^] Among AzTEC 1-15, only two sources (or 
13 ± 9%), AzTECI 1 and AzTEC14, were resolved into two dis¬ 
tinct components by Younger et al. (2007, 2009) in their 2" res¬ 
olution SMA 890 /mi imaging. We note that the northern and 
southern components of AzTEC 11 could just belong to an ex¬ 
tended object ( [Koprowski et al. 2014| . The 890 jjm flux den- 

8 In most cases, namely AzTEC17, 19, 21, 24, and 26, most of the 
components exhibit comparable 1.3 mm flux densities, i.e. the single¬ 
dish measured flux density from these sources appears to include con¬ 
tributions from equally bright components. 
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Fig. 4. Spectral energy distributions of the optical-to-IR counterparts of our SMGs. The red filled circles with error bars represent the photometric 
data points, while upper limits are denoted by downward pointing arrows [some of the data points with large error bars, such as the GALEX NUV 
(T e ff - 0.2316 jam), are also upper limits]. The solid lines represent the best-fit HYPERZ model to these data from the spectral model library of 
Michalowski et al. (2010). The panels on the right side show the corresponding total x 2 distributions as a function of redshift. The number of 
degrees of freedom (dof) in the^ 2 minimisation is indicated in the top right comer of each panel. The grey shaded area in the AzTEC 18, 19a, 
21a, 21b, 23, and 26a panels indicates the z range ignored in the determination of the best-fit photo-z solution (two almost equally probable z pho t 
solutions, where the higher z pho t is supported by the radio non-detection; see Appendix C). The thick red vertical line marks the best-fit z pho t value, 

shows the 99% CL The vertical dashed line in the AzTEC17a panel marks the spectroscopic redshift z spec = 0.834 

(Appendix C). ’ I—I 
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Fig. 4. continued. 


sity ratio for the two components of AzTEC 11 is 44 ± 23%, 
and even higher, 77 ± 35%, for AzTEC 14 ( Young er et al. 2009 


Table 1 therein). The low observed multiplicity fraction among 
AzTEC 1-15 could be the result of a sensitivity too low to reveal 
the real multiplicity (see Wang et al. 201 1| . However, as our an¬ 
gular resolution is only slightly better (and the observing wave¬ 
lengths are quite similar, i.e. 1.3 mm compared to ~ 0.9 mm), our 
observations provide a hint that the multiplicity fraction could 
be somewhat higher among the fainter SMGs AzTEC 16-30 (i.e. 


~ 30 - 40% compared to 13%). Considering the combined sam¬ 
ple AzTEC 1-30, the multiplicity fraction, based on the currently 
available data, appears to be ~ 25 ± 9%. 


Besides the present work, it has been found that interfer¬ 
ometric observations have the potential to resolve SMGs into 
separate components (e.g.[Daddi et al. 2009a; Wang et al. 2011 
ISmolcic et al. 2012bl |Barger et al. 2012[ IKarim et al. 2013 


Hodge e t al. 2013| ). For example, Smolcic et al. (2012b) used 
PdBI 1.3 mm observations at ~ 175 resolution of 28 SMGs 
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Fig. 4. continued. 


in the COSMOS field in conjunction with earlier interfero- 
metrically identified COSMOS SMGs, and concluded that 
> 15%, and possibly up to ~ 40% of single-dish detected 
SMGs (at 18" with AzTEC and at 2776 with LABOCA), 
consist of multiple sources. Consistent with this, Hodge et al. 
(2013) found that 24 out of their sample of 69 SMGs (~ 35%) 
detected with LABOCA at 870 pm (19" resolution) are split into 
multiples when observed with ALMA at the same wavelength 
[the ALMA-identified SMGs from the LABOCA Extended 





Chandra Deep Lield South (ECDLS) Submillimetre Survey 
(LESS), i.e. the AL ESS SMGs; see also I Karim et al. 20l3l 
|Simpson et al. 2014| . We also exa mined the mul tiplicity frac¬ 
tion among those LESS SMGs ( [Weifi et al. 2009| ) that have 
LABOCA 870 pm flux densities corresponding to our AzTEC 
1.1 mm flux density range, i.e. 3.3 mJy < S^^ m m < 9.3 mJy. 
Assuming that /3 = 1.5, this flux density range is 7.7 mJy 
< SsiOum <21.8 mJy. Based on their main and supplementary 
source samples, altogether 20 SMGs from Hodge et al. (2013) 
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have LAB OCA 870 p m flux densities in this range. Among 
these SMGs, five were found to have multiple (2-3) compo¬ 
nents (ALESS SMGs), resulting in a multiplicity fraction of 
25 ± 11%, which is very similar to our value. Dividing these 
LESS SMGs into two subsamples corresponding to the flux 
densities of AzTECl-15 and 16-30 (i.e. six SMGs with 9.8 
mJy < S' 87 Ofim ^ 21.8 mJy, and nine SMGs with 7.7 mJy 
< S 870 jum <9.1 mJy), we derived the corresponding multiplicity 
fractions to be 50 ± 29% and 44 ± 22%. These two values are 
similar within the counting uncertainties and hence, unlike what 
we found among AzTECl-30, the fainter LESS SMGs do not 
appear to exhibit a higher multiplicity fraction than the brighter 
SMGs. 

As recently discussed by Koprowski et al. (2014), there is 
some controversy about how common SMG multiplicity actu¬ 
ally is. The multiplicity statistics reported in the literature so far 
seem to suggest that the fraction of single-dish detected SMGs 
being composed of more than one SMG can be quite high (val¬ 
ues ranging from ~ 15% to ~ 40%). The multiplicity fraction 
also depends on the angular resolution of both the single-dish 
observations of the initial SMG detection and the follow-up ob¬ 
servations (the higher the former is the lower the multiplicity 
fraction is expected to be for a given follow-up resolution). How¬ 
ever, whether it is a common feature (several tens of percent) is 
an important knowledge when studying the number counts of 
SMGs, and future high-resolution observations of large, well- 
defined samples of SMGs are required to better understand the 
multiplicity fraction of submm-emitting galaxies. The completed 
ALMA array is well suited for this purpose. 

5.2. Counterpart associations of the SMGs AzTEC16-30 

Some of the detected (candidate) SMGs appear to have no coun¬ 
terparts at optical-to-IR wavelengths. These include AzTEC20, 
22, 24a-c, 26b, 28, and 30 [moreover, the I -band sources ly¬ 
ing 0735 from AzTEC 16, 1715 from AzTEC27, and 1741 from 
AzTEC29a might be unrelated to the (candidate) SMGs be¬ 
cause no sources at other wavelengths are identified there]. 
In particular, the case of AzTEC28, a clearly detected PdBI 
1.3 mm source (S/N = 5.5) without multiwavelength coun¬ 
terparts, shows that SMGs can be so heavily obscured by dust 
and/or at high redshift that only FIR-to-mm continuum emis¬ 
sion can be detected (this is probably true also for AzTEC27). 


Given the deep multiwavelength data available for the COSMOS 
field, for example the 1st UltraVISTA data release (DR1) go¬ 
ing down to K s < 24 mag (IMcCracken et al. 2012L the fraction 
of sources that lack shorter-wavelength (and radio) counterparts 
seems high (8/22 ~ 36%, or 11/22 = 50% if the three additional 
SMGs having only a nearby ACS I -band source are counted). 
More sensitive observations could reveal the presence of faint 
emission at optical-to-MIR waveleng ths, suc h as_ UltraVISTA 
DR2 reaching K s ~ 25 mag (IMcCracken et al. 2013b . and new 
IRAC 3.6 and 4.5 ji m observations (reaching ~ 25.5 mag) 
from the Spitzer Large Area Survey with Hyper-Suprime-Cam 
(SPLASH) (PI: P. Capak; ISteinhardt et al. 20l"4b . Younger et al. 
(2009) found that AzTEC13 and AzTEC14-E and 14-W are not 
coincident with any optical, Spitzer , or VLA sources. Therefore, 
among AzTEC 1-15, altogether comprising 17 SMGs, the frac¬ 
tion of SMGs that lack multi wavelength counterparts is 3/17 
or 18%, which is lower than for the fainter SMGs AzTEC lb- 
30. For comparison, 45% of the ALESS SMGs were found 
to lack MIR/radio counterparts ( Hodge et al. 201 3[ see also 
|Simpson et al. 2014"} cf. |Biggs et al. 2011 ). This fraction is the 
same (~ 45%) if we consider the 20 LESS SMGs that have sim¬ 
ilar flux densities as our AzTEC SMGs (see Sect. 5.1); the to¬ 
tal number of ALMA-detected components in these SMGs is 22 
(main and supplementary samples in |Hodge et al. 2013| ), and 12 
of them were found to have robust MIR/radio counterparts. 

The lack of multiwavelength counterparts means that we 
are not able to place strong constraints on the source’s photo¬ 
metric redshift. In particular, the absence of an optical counter¬ 
part suggests that the source is highly obscured by interstellar 
dust, which conforms to the fact that SMGs represent very dusty 
galaxies. More generally, the non-detections at optical, NIR, and 
cm wavelengths suggest that the source lies at a high redshift 
(z > 3 - 4), so that the flux density at a wavelength other than 
(sub)mm dims because of a large luminosity distance (i.e. the 
radiation suffers from the positive /^-correction). In contrast, the 
(sub)mm flux density stays almost the same over the redshift 
range z ~ 1-8 because of the negative /^-correction of the ther¬ 
mal dust emission ( [Blain & Longair 1993| . 

As shown in Table [3j for AzTEC 17a, 17b, and 23, the 
projected angular offset between the optical-to-NIR candidate 
counterpart (from the COSMOS/UltraVISTA catalogue) and the 
1.3 mm emission peak is quite large, 1724 - 1738. Moreover, 
for AzTEC 18, 26a, and 29b the above offset is also relatively 
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Table 4. Best available redshifts for the 30 brightest JCMT/AzTEC- 
detected SMGs. 


Source 

Redshift 

Comment and 



reference 


AzTEC1 
AzTEC2 a 


AzTEC3 


AzTEC4 

AzTEC5 


AzTEC6 a 


AzTEC7 

AzTEC8 

AzTEC9 fl 


AzTEC 10° 


azteciP 

AzTEC 12 
AzTEC 13 


AzTEC 14-E C 


AzTEC 14-W c 
AzTEC15 


4.3415 

1.125 

(4.28 ± 0.82) 
(3-60!°;!D 

5.298 

4 Q^+0.43 

P 

o.uj_ 028 
(1.85 + 0.23) 
>3.52 

(3-86^-91) 

2.30 + 0.10 
3.179 

1 07 +on 
A,u/ -o.io 

(1.357) 
(2.82 ± 0.76) 
(4.85!°;S) 

2 79 + 1-86 


1.599 
9 C4+0.13 
^ -0.33 

>4.07 

(4-VOflS) 

> 2.95 


1 30 +012 


spec-z; Yun et al., in prep, 
spec-z; Balokovic et al., in prep, 
radio/submm-z; this work 
radio/su bm m-z; 
|Koprowski et al. 2014 


spec -z:IRiechers et al. 20 




[Capaket al. 2011 1 
photo-zfSmolcic et al. 2012b 
photo-z: ISmolcic et al. 2012bl 
radio/submm-z; this work 
radio/submm-z; this work 
radio/su bmm-z; 

|Koprowski et al. 2014} 
photo-z: ISmolcic et al. 2012bl 
spec-z; Riechers et al., in prep. 
photo-z: ISmolcic et al. 2012bl 
spec-z; Salvato et al., in prep, 
radio/submm-z; this work 
p hoto-z; _ 

|Koprowski et al. 2014 


photo-z: ISmolcic et al. 2012b1 


photo-z; 

|Koprowski et al. 2014| 
spec-z; Salvato et a l., i n prep . 
photo-z; ISmolcic et al. 2012b! 
radio/submm-z; this work 
radio/ submm-z; 

|Koprowski et al. 2014| 
radio/submm-z; this work 
radio/sub mm-z; 

|Koprowski et al. 2014| 
photo-z; ISmolcic et al. 2012b] 
photo-z;ISmolcic et al. 2012bl 


AzTEC16 

> 2.42 

radio/submm-z; this work 

AzTEC17a 

0.834 

spec-z; Salvato et al., in prep. 


(0-75tg) 

photo-z; this work 


(2.29 ± 0.42) 

radio/submm-z; this work 

AzTEC 17b 

4.14 + °- 87 

photo-z; this work 


(> 2.49) 

radio/submm-z; this work 

AzTEC18 

3.00!g;l? 

photo-z; this work 


(> 2.20) 

radio/submm-z; this work 

AzTEC19a 

3 20 +018 
J,z,u -0 45 

photo-z; this work 


(4.22 ± 0.91) 

radio/submm-z; this work 

AzTEC19b 

1.11+0.10 

photo-z; this work 


(> 6.57) 

radio/submm-z; this work 

AzTEC20 

>2.35 

radio/submm-z; this work 

AzTEC21a 

2.60!°;| 8 

photo-z; this work 


(> 3.45) 

radio/submm-z; this work 

AzTEC21b 

2-8022:1 % 

photo-z; this work 


(> 2.47) 

radio/submm-z; this work 

AzTEC21c 

> 1.93 

radio/submm-z; this work 

AzTEC22 

>3.00 

radio/submm-z; this work 

AzTEC23 

1 60 +0 - 28 

photo-z; this work 


(> 2.06) 

radio/submm-z; this work 

AzTEC24a 

> 2.35 

radio/submm-z; this work 

AzTEC24b 

> 2.28 

radio/submm-z; this work 

AzTEC24c 

>3.17 

radio/submm-z; this work 

AzTEC26a 

2 50 +0 - 24 
Z ' JU - 0 . 14 

photo-z; this work 


(> 1.87) 

radio/submm-z; this work 

AzTEC26b 

> 1.79 

radio/submm-z; this work 

AzTEC27 

>4.17 

radio/submm-z; this work 

AzTEC28 

> 3.11 

radio/submm-z; this work 

AzTEC29a 

>2.96 

radio/submm-z; this work 

AzTEC29b 

1 45 +0 - 79 
-0.38 

photo-z; this work 


(> 7.25) 

radio/submm-z; this work 

AzTEC30 

> 2.51 

radio/submm-z; this work 


Notes. When multiple values are given for the redshift, the one not en¬ 
closed in parentheses has been adopted in the present study. Among 
AzTEC 16-30 the source names highlighted in bold-face indicate de¬ 
tections with S/Ni. 3mm > 5.5. (a) See text for details about the differ¬ 
ence between our redshifts and those from Koprowski et al. (2014). 

AzTEC 11 was resolved into two 890 pm sources (N and S) by 
Younger et al. (2009). Here, for the redshift analysis, we treat it as 
a single source becaus e the two components are probably physically 
related ( [Koprowski et al. 2014]) . (c) AzTEC 14 was resolved into two 
890 pm sources (E and W) by Younger et al. (2009). The eastern 
component appears to lie at a higher redshift than the western one 
dSmolcic et al. 2012bl . 
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Fig. 5. Top: The distribution of redshifts of the studied SMGs. The blue 
filled histogram shows the redshift distribution of the SMGs AzTEC 1- 
15, while the red filled histogram shows that of AzTEC 16-30. The open 
grey histogram indicates the redshift distribution of the combined sam¬ 
ple (AzTEC 1-30). The vertical dashed lines mark the corresponding 
median redshifts (blue: z = 3.05; red: z = 3.20; grey: z = 3.17). The 
lower redshift limits (see Table [4j were placed in the bins correspond¬ 
ing to those values, but the indicated median redshifts were properly 
derived through survival analysis. Bottom: The probability density dis¬ 
tribution of the redshifts of AzTEC 1-30. The light grey curve shows the 
unsmoothed distribution, and the black curve represents the Gaussian- 
smoothed kernel density estimate (see text for details). The median red¬ 
shift and the 68% Cl are indicated. 


large (0782, 0794, and 0776, respectively). Such an offset could 
be the effect of complex source morphology, expected in the 
case of galaxy mergers ( Daddi et al. 200^, and/or be the resul t 
of strong differential dust obscuration (e.g. ICarilli et al. 20101) . 
With respect to our sample of 22 detections, ~ 14-27% (three 
to six sources) exhibit a considerable (> 078) projected sepa¬ 
ration between the PdBI and UltraVISTA emission peaks. Even 
though the positional error of our PdBI sources is much lower 
(~ 072; Sect. 2.2), the above mentioned angular separations 
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are still within the beam FWHM of ~ IV8. Moreover, towards 
AzTEC 17a, showing the largest offset between the PdBI peak 
emission and the UltraVISTA source in our sample (1V38), the 
Spitzer/ IRAC counterpart lies only 073 from the PdBI peak, 
strongly indicating that the SMG is emitting at observing-frame 
wavelengths > 3.6 pm. The photo-z value of 0.75 +q 22 we de¬ 
rived for AzTEC 17a is in good agreement with the spec-z of 
0.834, further strengthening our counterpart identification (Ap¬ 
pendix C). 

5.3. Redshift distribution of the JCMT/AzTEC-detected 

SMGs in the COSMOS field, and comparison with other 
surveys 

The median redshift of the SMGs studied here (AzTEC 16-30) 
is found to be z = 3.20 ± 0.25, while that for the SMGs 
AzTECl-15 is z = 3.05 ± 0.44. The combined sample of these 
JCMT/AzTEC-detected SMGs, i.e. AzTEC 1-30, has a median 
redshift of z = 3.17 ± 0.27, which corresponds to an age of the 
universe of 2.06 +q 21 Gyr or about 15+ 2 % of its current age. A 
two-sample K-S test of the null hypothesis that the two subsam¬ 
ples, AzTEC 1-15 and AzTEC 16-30, are drawn from the same 
underlying parent distribution yielded a p-value of 0.7342. This 
suggests that the two subsamples are probably sampled from a 
common distribution. The /-test also suggests that the mean red- 
shifts of the subsamples are similar to each other. In Fig. [6} we 
show the redshift distribution of AzTEC 1-30 and, for compari¬ 
son, those derived in other SMG surveys discussed below. 

We note that our combined SMG sample contains a source 
(AzTEC 17a) at a redshift of ^ 0.8. This redshift is quite similar 
to that of the ~ 25 Mpc long filamentary COSMOS large-scale 
structure (the COSMOS Wall) at z - 0.73 ( Guzzo et al. 2007) . 
However, cross-correlation with the redshift survey catalogue 
consisting of 1 023 galaxies belonging to the COSMOS Wall 
did not yield any cross matches within a 175 search radius 
(A. Iovino, priv. comm.). Although AzTEC17a appears to be a 
member of a z ~ 0.8 galaxy overdensity (V. Smolcic et al., in 
prep.), our redshift survey is not subject to strong cosmic vari¬ 
ance arising from the COSMOS large-scale structure, and the 
different results compared to other cosmological survey fields 
have their origin in other effects (e.g. observing wavelength used, 
inclusion of radio-faint SMGs, etc.). 

As demonstrated in the present study, radio-dim SMGs 
are probably lying at high redshifts (cf. |Chapman et al. 2005) . 
For example, the z - 4.3 SMG AzTEC 1 is associated (near 
the SMA position) with only a weak 20 cm radio source 
(*5 20 cm = 48+14 pjy), and th e z - 5.3 SMG AzTEC3 has 
no 20 cm counterpart ( [Younger et al. 2007| . These results are 
based on the VLA 1.4 GHz imaging down to a mean lcr rms 
depth of - 10.5 pjy beam -1 (ISchinnerer et al. 20071 . However, 
both AzTEC 1 and AzTEC3 are associated with VLA 10 cm 
emission where the corresponding maps have a lcr noise of 
4.5 pjy beam -1 , which, for a typical radio spectral index of 
-0.7, corresponds to the equivalent 20 cm rms noise of about 
~1.4 times lower than the aforementioned 1.4 GHz sensitivity 
level (ISmolcic et al. 20151) . For comparison, the 1 cr rms noise at 
1.4 GHz in the phase centre of the seven fields analysed by Chap¬ 
man et al. (2005) was 4-15 ply beam -1 , while the rms sensitivity 
in the 1.4 GHz imaging of the ALESS SMGs was 6 pjy beam -1 
( Thomson et al. 2014 ). 

Some earlier studies of SMGs have suggested that more lu¬ 
minous SMGs lie, on average, at higher redshifts compared to 
less luminous SMG s (e.g. Ilvison et al. 20021 Pope et al. 2 006; 
|Younger et al. 2007] |Biggs et al. 2011| ISmolcic et al. 2012bi 


Our redshift analysis suggests that the brighter SMGs (AzTEC 1- 
15) have a similar mean redshift ((z) = 3.16 ± 0.37) within the 
errors to the less bright sources in our sample (AzTEC 16-30; 
(z) = 3.02 ± 0.20). Furthermore, as noted above, the median 
redshifts of the two subsamples are similar to each other within 
the uncertainties. Wall et al. (2008) suggested that there might 
be two SMG subpopulations, divided by their luminosity (divid¬ 
ing line being at LssOyum = 3 x 10 12 L 0 ): these would evolve in 
different ways, and the corresponding luminosity functions have 
different shapes. However, in agreement with our result, Ward- 
low et al. (2011) found no significant correlation between the 
redshift and 870 pm submm flux density for their extensive sam¬ 
ple of LESS SMGs, although it should be noted that many (~ 
35% - 50%) of the LES S SMGs have been resolved into multi - 
ple sources with ALMA (IKarim et al. 20131 Hodge et al. 2013 ). 
Moreover, as discussed by Hodge et al. (2013; Sect. 5.5 therein), 
some of the LESS SMGs suffered from missed/misidentified 
multi wavelength counterparts, which means that they had incor¬ 
rect photometric redshifts. Overall, - 45% of the ALESS SMGs 
were missed by the sophisticated counterpart association util¬ 
ising multiwavelength information by Biggs et al. (2011), and 
of the reported counterparts - 1/3 were found to be incorrect 
( [Hodge et al. 2013[ ). However, the ALESS SMGs also do not 
exhibit any significant trend between the redshift and S 870 yum 
( [Simpson et al. 2014] their Fig. 7). The recent semi-analytic 
modelling of 850 pm SMG surveys by Cowley et al. (2015) pre¬ 
dicted that the bright SMG population (S 850 yum > 5 mJy) has a 
lower median redshift than the faint SMG population (1 mJy < 
S 850 yum < 5 mJy). We note that all our 1.1 mm SMGs AzTEC 1- 
30 would belong to the bright SMG population of Cowley et al. 
(2015), i.e. all our SMGs have Si.imm > 2 mJy when scaling 
the S'85Oyum > 5 mJy threshold by assuming that = 1.5. From 
the 50 mock surveys of 0.5 deg 2 in angular size by Cowley et 
al. (2015), where SMGs were generated out to z = 8.5, the me¬ 
dian redshift for the bright SMGs was derived to be 2.05, while 
that for the faint SMGs was found to be 2.77. The authors also 
compared their model predictions with the ALESS SMG survey, 
and found that the model successfully reproduces the median 
redshift of the ALESS photo-z distribution (see below). The op¬ 
posite redshift trend predicted by Cowley et al. (2015) compared 
to some earlier observational results can, as suggested by the au¬ 
thors, be tested with future interferometric SMG surveys. They 
also pointed out that field-to-field variance can play a role when 
comparing theoretical model predictions with observational sur¬ 
vey results. 

In the following, we investigate the origin of differences 
in mean/median redshift for differently selected SMG samples. 
WeiB et al. (2013) carried out a blind redshift survey with 
ALMA towards 26 strongly lensed SMGs originally detected 
with the 10 m South Pole Telescope (SPT) at 1.4 mm. Their 
sample consisted of sources with high 1.4 mm flux densities of 
> 20 mJy, and the average redshift of the sample was found 
to be (z) = 3.5, significantly higher than what is found for 
radio-identified SMGs, but quite similar to that of AzTEC 1- 
30 selected at 1.1 mm ((z) = 3.19 ± 0.22). Simpson et al. 
(2014) presented the first photo-z distribution for the ALESS 
SMGs derived using HYPERZ SED fitting with the spectral tem¬ 
plates of Bruzual & Chariot (2003) rather than those optimised 
for SMGs by Michalowski et al. (2010) we have used (see 
ISmolcic et al. 2012al b for the comparison of these model li¬ 
braries in the photo-z analysis). For their sample of 77 SMGs 
with broadband photometry, Simpson et al. (2014) found the 
median redshift to be z = 2.3 ±0.1 (2.5 + 0.2 when the 19 
sources with poorer photometry were included in the analy- 
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sis). This is very similar to the median spectroscopic redshift 
of z = 2.2 derived by Chapman et al. (2005) for a sample of 
73 radio-identified (VLA 1.4 GHz) SCUBA SMGs compiled 
from seven separate fields, but lower than the median redshift 
z = 3.17 ± 0.27 we derived for AzTECl-30. To perform a more 
quantitative comparison with the ALESS SMG redshift distri- 
bution, we selecte d only those LESS SMGs ( [Weifi et al. 2009[ 
|Hodge et al. 2013| ) with LAB OCA 870 pm flux densities cor¬ 
responding to our AzTEC 1.1 mm flux density range (see 
Sect. 5.1). Altogether 17 ALESS SMGs from Simpson et al. 
(2014) satisfy this criterion (7.7 mJy < 5 870 yum < 21.8 mJy). For 
this 870 pm flux limited sample we derived the following photo¬ 
metric redshift statistics: (z) = 2.99 ± 0.31, z = 2.85 ± 0.39, 
std = 1.27, and Cl = 2.39 - 3.59 (95%). As can be seen the 
median value is higher than that for the original sample of 77 
ALESS SMGs ( [Simpson et al. 2014] ), but lower than the me¬ 
dian redshift for AzTEC 1-30. We also carried out a K-S test 
between our sample (excluding the lower limits in which case 
(z) = 2.61 ± 0.26) and the ALESS SMGs with comparable flux 
densities, and found that they could have a common underlying 
parent distribution (D = 0.2379, p = 0.6379). Furthermore, the 
Mest suggests that these two samples have a comparable aver¬ 
age redshift (p = 0.3573 under the null hypothesis that the two 
(z) values are equal). The facts that the LESS and JCMT/AzTEC 
SMGs were selected at different wavelengths (870 pm compared 
to 1.1 mm) from different fields, and that the ALESS sample 
is larger than ours make a direct comparison between the two 
questionable. Given that the lower redshift limits to some of our 
SMGs raise the total sample average to (z) = 3.19 ± 0.22 could 
be an indication that the 1.1 mm wavelength selects somewhat 
higher-redshift SMGs. Moreover, the ALESS sample was drawn 
from the Extended Chandra Deep Field South that has a size of 
30' x 30', or 0.25 m°, while our SMGs were selected from the 
0.15 m° COSMOS JCMT/AzTEC field. The intrinsic target field 
properties, or cosmic variance might therefore also play a role 
(cf. Cowley et al. 2015| ). 


Swinbank et al. (2014) analysed the ALESS SMGs detected 
in at least two Herschel /SPIRE bands. They found that the 
SMGs exhibiting the peak dust emission at A = 250, 350, and 
500 pm have median redshifts of z = 2.3 ± 0.2, 2.5 ± 0.3, and 
3.5 ± 0.5, respectively (the SPIRE non-detected sources were 
found to have a median photo-z of 3.3 ± 0.5). Although the SED 
peak position is not always well constrained_and the s ubsa mples 
clearly overlap with each other (Fig. 6 in lSwinbank et al. 2014k 
this suggests that there is a positive correlation between the 
source redshift and the SED peak wavelength. Within the er¬ 
rors the 500 pm peakers have a median redshift comparable to 
that of our SMGs. The highest redshift SMG known to date, 
HFLS3 at z - 6.34, is also a 500 pm peaker and was orig¬ 
inally found from the HerMES survey as having a very high 
500 yum flux density of 5 500 yum = 1.46x5 350 ^ = 47.3 + 2.8 mJy 
dRiechers et al. 20131) . Similarly, Dowell et al. (2014) selected 
dusty star-forming galaxies from the HerMES survey on the 
basis of their Herschel /SPIRE colours, and found most of the 
500 pm peakers (including HFLS3) to lie at very high redshifts 
(z > 4). The relatively high redshifts among the AzTEC 1-30 
COSMOS SMGs could be a selection effect in the sense that they 
were originally discovered at A = 1.1 mm, although cosmic vari¬ 
ance can also play a role because the COSMOS field is known to 
contain a relatively large number of very high-z SMGs. Zavala et 
al. (2014) carried out simulations of the SMGs’ redshift distribu¬ 
tions, and they studied how different selection effects affect the 
derived distributions. Their simulated data showed the increase 
in the median redshift as a function of wavelength (changing 
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from z = 2.06 ± 0.10 at 450 pm to z = 2.91 ± 0.12 at 2 mm). 
However, they demonstrated that the differences reported in the 
literature can be explained by the observing wavelength (related 
to the SED temperature) used and, to a lesser degree, by the map 
noise level and angular resolution, and that some of the redshift 
distributions suggested to be different from each other can actu¬ 
ally be drawn from the same parent distribution. 

As discussed above, the derived SMG redshift distribution 
can be highly affected by the wavelength selection and source 
sample under study. A well established example is the radio pre¬ 
selection that biases the sample towards lower-redshift (z < 3) 
SMGs ( [Chapman et al. 2005] ). However, spectroscopic observa¬ 
tions have shown that the z > 4 SMGs are more common than 
originally thought (see references in Sect. 1). A mix of differ¬ 
ent methods to derive the redshifts, such as spectroscopic and 
photometric method can also lead to a biased distribution of red¬ 
shift values. For galaxies in the COSMOS field, however, it has 
been shown that the photo-z values agree well with those derived 
through spectroscopic observations (I lbert et al. 2013 ). Consid¬ 
ering only the most secure spec-z values at K s < 24 (a sample 
of 12482 galaxies), Ilbert et al. (2013) found that the photo- 
z accuracy is <ta z /(i+ z ) = 9.6 x 10 -3 and only 2.1% are catas¬ 
trophic failures with |z p hot _ £specl/(l + z spe c) >0.15. The different 
methods of deriving the photo-z values (e.g. varied assumptions 
and spectral templates) can also lead to differing results, but our 
photo-z values derived from HYPERZ using the SMG SED tem¬ 
plates from Michalowski et al. (2010) are expected to be reliable 
(ISmolcic et al. 2012al b): among our new SMG sample, this is 
supported by the case of AzTEC 17a (z spe c = 0.834 compared 
to zphot = 0-25 12 )- I n some cases the best photo-z solution 
is uncertain because the corresponding xlot distribution is com¬ 
plex having a broad minimum or multiple dips of comparably 
low^ ot value. Moreover, our SMG redshift distribution is partly 
based on lower limits only that were derived using the Carilli- 
Yun radio-submm redshift indicator dCarilli & Yun 19991 2000). 
This method is subject to a degeneracy between T d us t and z, and 
can suffer from large uncertainties. Another caveat in determin¬ 
ing the photo-z values is the possible contamination by AGNs. 
The reason for this is that methods based on stellar libraries 
might suffer from short-wavelength (UV to MIR) AGN emis¬ 
sion (see Smolcic et al. 2012bl for further discussion). However, 
as mentioned earlier our sources do not exhibit any strong X-ray 
signatures and are therefore unlikely to contain bright AGNs. 

To summarise, our new interferometric observations have en¬ 
abled us to pinpoint the multiwavelength counterparts of our 
SMGs, and therefore to derive the photo-z values for these 
SMGs. For this type of analysis, interferometry provides an im¬ 
portant improvement because the usage of single-dish (sub)mm 
data of ~ 10 - 30" resolution can result in a wrong counterpart 
identification, and therefore also wrong redshift of the SMG. For 
five SMGs among AzTEC 1-15, we have a secure spectroscopic 
redshift available, but only one spec-z among AzTEC 16-30. In 
the ideal case, all the SMG redshifts would be based on spec¬ 
troscopic data. This way one could carry out a completely fair 
comparison between our two subsamples of AzTEC 1-15 and 
AzTEC 16-30. 


6. Summary and conclusions 

We have used the IRAM/PdBI to carry out an interferometric 
1.3 mm continuum follow-up study of a sample of 15 SMGs 
originally detected in the COSMOS field with the JCMT/AzTEC 
bolometer at 1.1 mm (51.1mm - 3.3 - 3.9 mJy) by Scott et al. 
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Fig. 6. The SMG redshift distributions discussed in the text. Besides our 
result for AzTEC 1-30, the plot shows the z distributions for the radio- 
identified SCUBA SMGs from Chapman et al. (2005), ALESS SMGs 
from Simpson et al. (2014), and SPT SMGs from WeiB et al. (2013). The 
vertical dashed lines show the corresponding median redshifts [z = 3.17 
for AzTEC 1-30, z- 2.2 for the Chapman et al. (2005) SMGs, z = 2.3 
for the ALESS SMGs, and z = 3.7 for the SPT SMGs]. 

(2008). The good angular resolution of about 178 allowed us 
to pinpoint the positions of the actual SMGs giving rise to the 
single-dish detected mm emission. We combined these new ob¬ 
servations with results from the literature to study the ensemble 
properties of the 30 most significant (S/N = 4.0- 8.3) SMGs de¬ 
tected by JCMT/AzTEC in the COSMOS field. Our main results 
and conclusions are summarised as follows: 

1. The total number of sources detected in this survey is 22, 
where the sample does consist of S/Ni. 3 mm > 4.5 detections 
(whether or not having a counterpart) and sources detected 
with 4 < S/N 1<3mm < 4.5 that have multiwavelength coun¬ 
terparts. AzTEC 19 is found to be the most significant 1.3 mm 
emitter (in the observed frame) with S/Ni. 3 mm = 10.3. 

2. Visual inspection of the 1.3 mm images reveals that 
AzTEC19, 21, 27, and 28 have elongated/clumpy morpholo¬ 
gies, a possible manifestation of galaxy merging. AzTEC27 
appears to be a gravitationally lensed SMG, where two in- 
tervenin g gala xies are warping and magnifying the radiation 
(see Fig. & 

3. AzTEC17, 19, 21, 24, 26, and 29 appear to split into two 
to three sources of 1.3 mm emission. Considering our new 
SMG sample (15 single-dish detected sources), this would 
mean that the multiplicity fraction is 40 ± 16% (33 ± 
15% if AzTEC24 is spurious). In all the other cases ex¬ 
cept AzTEC29, the source components have comparable 
observing-frame 1.3 mm flux densities. Among AzTEC 1- 
15 there are two SMGs that are known to be multiple sys¬ 
tems. Combining these statistics we conclude that the mul¬ 
tiplicity fraction among AzTECl-30 is ~ 25 ± 9%. Deep, 
high-resolution (sub)mm surveys of large SMG samples are 
required to unambiguously determine how common multi¬ 
plicity is among SMGs. 

4. Besides the spectroscopic redshift of AzTEC 17a, the red- 
shifts of AzTEC 16-30 were derived using either opti- 
cal/IR photometric data or the Carilli-Yun redshift indicator 


( Carilli & Yun 19991 , 2000). In many cases only lower lim¬ 
its could be estimated, and the median redshift was found to 
be z = 3.20 + 0.25. We identified some high-redshift can¬ 
didates; in particular, AzTEC 17b has a photo-z of 4.14^^, 
and a lower limit to Zradio/submm of AzTEC27 was derived to 
be as high as > 4.17. For the 15 brightest JCMT/AzTEC 
1.1 mm detected SMGs, namely AzTEC 1-15, the median 
redshift is z = 3.05 ± 0.44 (partly based on secure spectro¬ 
scopic redshifts). For the combined sample of AzTEC 1-30, 
the median redshift was found to be z = 3.17 ± 0.27. This is 
higher than what is usually reported for SMGs, but in agree¬ 
ment with mm-selected SMG samples. 

5. The redshift analysis does not support the earlier obser¬ 
vational result that brighter SMGs (our sources AzTEC 1- 
15) would lie at higher redshifts than the fainter SMGs 
(AzTEC 16-30). Instead, besides the median redshifts, the 
mean redshifts of AzTEC 1-15 and AzTEC 16-30 are simi¬ 
lar to each other within the errors ((z) = 3.16 + 0.37 and 
(z) = 3.02 + 0.20, respectively). The /-test also supports the 
similarity between the mean redshift values. Finally, the K- 
S test suggests that the two subgroups are probably drawn 
from a common parent population, but we note that the high¬ 
est redshift (z > 4.3) SMGs are found among the strongest 
millimetre emitters. The absence of any significant trend be¬ 
tween the source redshift and millimetre flux density is in 
agreement with that_found for the ALESS SMGs at 870 jm l 
wavelength ( [Simpson et al. 2014| . 

Some of the great challenges in detailed observational stud¬ 
ies of SMGs is to reliably identify their multiwavelength coun¬ 
terparts. While interferometric (sub)millimetre imaging is a pre¬ 
requisite for secure counterpart identifications, faint SMGs, flux- 
boosted in single-dish observations, might not be detected in 
shallow interferometric maps. Moreover, as the present study 
demonstrates, the sample might suffer from spurious sources that 
may or may not have nearby sources detected at other wave¬ 
lengths. The knowledge of secure counterparts is required to 
obtain accurate estimates of the sources’ photometric redshifts. 
Spectral line observations of atoms (such as the T res t = 158 /i m 
forbidden C + fine-structure line) or molecules (particularly rota¬ 
tional lines of 12 CO) are needed to obtain the most reliable red¬ 
shifts [cf. ALMA observations towards SPT SMGs by WeiB et al. 
(2013)], and hence to push forward our understanding of high- 
redshift, submillimetre-luminous starburst galaxies, their role in 
galaxy formation and evolution, and cosmic star formation his¬ 
tory. 
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Appendix A: Multiwavelength images 

A selection of zoomed-in multiwavelength views towards 
AzTEC 16-30 is shown in Fig. |A.l| In the first (top left) panel of 
each source we show the PdBI 1.3 mm image overlaid with the 
same contour levels as in Fig.[l] The PdBI images are also anno¬ 
tated with the source designations. The positive 1.3 mm contours 
are overlaid on the other wavelength images to guide the eye. 

Appendix B: Updated redshifts of the 15 brightest 
JCMT/AzTEC-detected SMGs: AzTEC1-15 

Here we provide the reader with an overview of the redshifts 
of the SMA-detected SMGs AzTECl-15. Among these SMGs, 
there are eight spectroscopic redshifts reported in the literature: 
for AzTECl, 2, 3, 5, 6 , 8 , 9, and 11 (see ISmolcic et al. 2012bt 
their Tables 1 and 4 and references therein). Despite the partially 
extensive efforts and data coverage some of these redshift de¬ 
terminations are still uncertain. We discuss below the updated 
redshifts among AzTEC 1-15 and the cases where there is some 
confusion about the source redshifts in the literature. 

Smolcic et al. (2011) determined a spectroscopic redshift of 
4.650 ± 0.005 for AzTECl. The UV-NIR photometric redshift 
they derived, z p hot = 4.64^ ^, was found to be very similar to the 
Zspec value, although a secondary photo-z solution at z p hot = 4.44 
was also found. A somewhat lower photo-z of 4.26 +q2q was 
derived by Smolcic et al. (2012b) using the same method as 
in the present paper (Sect. 4.2). The CO spectral-line observa¬ 
tions using the Redshift Search Receiver (RSR) on the Large 
Millimetre Telescope (LMT) performed by M. S. Yun et al. (in 
prep.) yielded a spec-z value of 4.3421 for AzTECl. Their SMA 
follow-up observations of C + emission yielded a line detection 
at Zspec = 4.3415, in very good agreement with the CO observa¬ 
tions. Since it is based on interferometric observations, this last 
redshift is adopted in the present work. We note that the new 
spec-z of AzTECl explains the non-detection of the CO(5 - 4) 
line emission by Smolcic et al. (2011) because their PdBI and 
Combined Array for Research in Millimetre-wave Astronomy 
(CARMA) observations covered the redshift ranges 4.56-4.76 
and 4.94-5.02. 

The optical spectrum observed with the Deep Extragalac- 
tic Imaging Multi-Object Spectrograph (DEIMOS) on the 10 m 
Keck II telescope towards AzTEC2 exhibits an emission feature 
that can be assigned to the [On] T3727 forbidden-line doublet 
at Zspec = 1.124, and the J - 2 - 1 rotational line of CO de¬ 
tected with CARMA suggests a similar redshift (z spe c = 1.126; 
M. Balokovic et al., in prep.). Following Smolcic et al. (2012b), 
we adopt the value z spe c = 1.125 as the redshift of AzTEC2. 
Koprowski et al. (2014) claimed that the target position of these 
spectral line observations was 174 away from the SMA 890 p m 
position ( [Younger et al. 2007] ). They concluded that the SMG 
lies at a redshift of 3.60^]^ derived from the radio/submm flux- 
density ratio because the radio source is only 074 away from the 
SMA position. The redshift we derived from the radio/submm 
flux-density ratio is z = 4.28+0.82. The latter difference emerges 
because Koprowski et al. (2014) based their calculation on the 
average z - 2-3 SMG spectral template derived by Michalowski 
et al. (2010), while we utilised the Carilli-Yun redshift indicator 
dCarilli & Yun 19991 2000) as described in Sect. 4.2. However, 
as shown in Fig. |E.1[ the Keck/DEIMOS slit was centred 0798 
from the SMA position, and the spectrum was extracted from the 
SMA peak of the SMG. Moreover, the EW-oriented DEIMOS 
slit did not cover the optically visible foreground galaxy on the 
southern side of AzTEC2 (UltraVISTA ID 232116, z ph ot = 0.34; 


cf. Fig. B1 in Koprowski et al. 2014| ). This implies that the red¬ 
shift of AzTEC2 is close to unity instead of the higher value 
proposed by Koprowski et al. (2014). 

The spec-z of AzTEC5 was previously reported to be z spe c = 
3.971 (ISmolcic et al. 2012bt . The Keck/DEIMOS slit position 


E.l| We note that the slit 
axies other than AzTEC5, 


and orientation are shown in Fig. 
does not include emission from gal; 
but the corresponding DEIMOS spectrum is of poor quality. 
Therefore, we adopt the photo-z of 3.05^28 f rom Smolcic et 
al. (2012b). For comparison, Koprowski et al. (2014) derived 
a photo-z of 4.19 ^q. The redshift we derived from the ra¬ 
dio/submm flux-density ratio using the Carilli & Yun (2000) for¬ 
mula (see Sect 4.2), 1.85 ± 0.23, is also lower than the value 
2 . 9 O+ 0 J 5 calculated by Koprowski et al. (2014). 

Koprowski et al. (2014) argued that the spectroscopic red¬ 
shift of AzTEC 6 , Zspec = 0.802, is uncertain because it is mea¬ 
sured to wards an optically visible object about 1" from the SMA 
position ( [Younger et al. 2007 1 , and that the submm/radio flux ra¬ 
tio of AzTEC 6 is inconsistent with a low redshift (they derived a 
value of z = 3.86^^ from the radio/submm flux-density ra¬ 
tio, while we derive the value z > 3.52 because AzTEC 6 is 
not detected at 20 cm). The photo-z value derived by Smolcic 
et al. ( 2012 b), z p hot = 0.82^-}^, is similar to the z spe c value, 
while Koprowski et al. (2014) reported a value of z p hot = 1-12 
for this object. The above-mentioned optically visible galaxy 
lies only 0766 from the SMA position, and 0712 away from the 
Keck/DEIMOS slit centre (see Fig. ED- The above spec-z value 
was derived from a high quality spectrum (flag 4; J. S. Kartaltepe 
et al., in prep.) extracted from a position that lies 0762 from the 
SMA position, and coincides with the optical galaxy. This im¬ 
plies that the spectral line emission originates in this foreground 
object as suggested by Koprowski et al. (2014). A redshift of 
Zspec = 0.802 indeed conflicts with the non-detection of AzTEC 6 
at 20 cm, and we therefore adopt the redshift z > 3.52. 

The DEIMOS spec-z of AzTEC9 was reported to be 
1.357, and its c orresponding photo-z was found to be 1.07^ }^ 
(ISmolcic et al. 2012bt . However, the spec-z value is based on 
a relatively weak spectrum (M. Salvato et al., in prep.), and 
is therefore quite uncertain. Koprowski et al. (2014) reported 
that the above redshift values refer to an object as far as about 
278 from the SMA position ( [Younger et al. 2009] ). Again, this 
is not the case, but the DEIMOS spectrum was extracted from 
the SMA position (the slit centre was offset from the SMA peak 
by 0745; see Fig. ET). Koprowski et al. (2014) stated that the 
submm/radio flux-density ratio of AzTEC9 is inconsistent with 
a redshift value close to unity [they derived z = 4.60^ ^, while 
our result based on the Carilli & Yun (2000) redshift formula is 
z = 2.82 ± 0.76]. They also derived a high photo-z of 4.85^ ^ 
for AzTEC9 (counterpart lying 0777 from the SMA position). 
There is a Spitzer /IRAC source 0762 from the SMA position, 
and the Wardlow et al. (2011) redshift formula gives a redshift 
of z - 2.75, which is similar to the redshift we inferred from the 
radio/submm flux-density ratio, but considerably lower than the 
redshifts derived by Koprowski et al. (2014). In the present study, 
we adopt the photo-z of 1.07^-JJ from Smolcic et al. (2012b) 
because the corresponding;^ distribution exhibi ts a clear min¬ 
imum at that value (see Fig. 6 in ISmolcic et al. 2012bl) . 

For AzTECIO, the photo-z derived by Smolcic et al. (2012b) 
is 2.79^ 29 , while Koprowski et al. (2014) determined a photo- 
z of 5.00^-^ for the optical/NIR source about 175 from the 
SMA position. The most up-to-date COSMOS spectroscopic- 
redshift catalogue gives a likely (quality flag 2) DEIMOS red¬ 
shift value of Zspec = 0.547 towards AzTECIO (only 07018 off- 
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Fig. A.l. Multiwavelength views towards AzTEC 16-30. The panels from top left to bottom right for each source are as follows: PdBI 1.3 mm, 
YLA 20 cm, VLA 10 cm, Spitzer 24 pm, Spitzer 8 pm, Spitzer 3.6 pm, UltraVISTA YJHK S colour composite, and HST/ACS I- band. The overlaid 
1.3 mm contours are as in Fig. |T| and positive 1.3 mm contours are shown in all panels. The synthesised beam of the PdBI data is shown in the 
bottom left corner in the first panel for each source. A scale bar indicating the l" projected length is shown in the PdBI panel, and the corresponding 
proper length [kpc] at the indicated redshift is also denoted (except when only a lower limit to z could be derived). The catalogue positions of the 
Herschel/SPIRE 250 pm sources are marked with plus signs in the PdBI images towards AzTEC 19, 20, and 24. The diamond symbol in the PdBI 
image towards AzTEC24 indicates the position of the ASTE/AzTEC 1.1 mm source AzTEC/C48 from Aretxaga et al. (2011). 


set from the SMA position; M. Salvato et al., in prep.). As illus¬ 
trated in Fig. EH the DEIMOS slit however picked up emission 
from a foreground galaxy at z p h 0 t - 0.51 (ID 302846 in the new 
UltraVISTA catalogue) that lies 0'.'97 NW of AzTEC 10. Since 
AzTEC 10 is not detected at 1.4 GHz, it appears to lie at a high 


redshift. In the present study, we adopt the photo-z from Smolcic 
et al. (2012b), but note that because of the multiple nearby coun¬ 
terparts of this source (three within 2") it is difficult to obtain 
accurate photometry for AzTEC 10. 
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Fig. A.l. continued. 


Both AzTEC 13 and AzTEC 14-E have neither optical nor 
IRAC counterparts, and we derived lower limits of z > 4.07 
and z > 2.95 for their radio/submm flux-ratio based redshifts 
[these differ from the values z > 3.59 and z > 3.03 derived by 
Smolcic et al. (2012b) because of the different assumptions that 
we used here]. These lower limits are consistent with the cor¬ 
responding values of Koprowski et al. (2014), i.e. z = 4.70+J'j^ 
and z = 3 . 38 + 054 , respectively. The updated COSMOS spec- 
z catalogue gives a high quality (flag 4) DEIMOS redshift of 
Zspec = 0.471 for a target that is 07015 away from the SMA po¬ 
sition of AzTEC 13 (M. Salvato et al., in prep.), where the spec¬ 


troscopic slit centre was positioned 0760 away from the SMA 
peak). This redshift is much lower than the other estimates men¬ 
tioned above. However, as shown in Fig. |E.l| there are two fore¬ 
ground galaxies lying 2722 SE and 2745 SW from AzTEC 13 
(UltraVISTA IDs 268116 and 268129 with the photo-z values of 
0.49 and 0.45, respectively); these could have contaminated the 
spectral line measurements, although they do not lie within the 
slit boundaries. A low redshift of AzTEC 13 would indeed be in¬ 
consistent with the radio non-detection (cf. AzTEC 6 ). For both 
AzTEC 13 and AzTEC 14-E, we adopt the redshifts derived from 
the radio/submm flux ratio (z > 4.07 and z > 2.95). 
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Fig. A.l. continued. 


Appendix C: Multiwavelength counterparts and 
redshifts of the SMGs AzTEC16-30 

Below we describe the multi wavelength appearances of our 
PdBI SMGs and provide notes of their redshifts. 

AzTEC16. The 5<x point-like 1.3 mm emission feature near 
the edge of the PB FWHM appears to have no counterparts at 
other wavelengths. Altogether four negative sources were found 
in this field with |S/N| = 4.0 - 6.2, three of which lie out¬ 
side the PB. Hence, AzTEC 16 could be spurious despite its rel¬ 
atively high significance. We note that 4'.'3 west of the target 
field centre, there is a ~ 3<x (35.6 pJy beam -1 ) VLA 20 cm 
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source, which also appears to be detected at Spitzer /IRAC and 
24 pm wavelengths. This source can be identified as the galaxy 
J095950.03+024416.5 from the COSMOS optical/NIR cata¬ 
logue ( [Capak et al. 2007| ); its UltraVISTA DR1 photometric - 
redshift catalogue ID is 319194 (z - 1.62; filbert et al. 20l3l ) . 
The radio non-detection of AzTEC 16 suggests a lower limit 
of z > 2.42 to its redshift. In Table 1 of Smolcic et al. 
(2012b), the source called AzTEC 16 at a spectroscopic red- 
shift of 1.505 (based on high-resolution CO observations with 
CARMA; K. Sheth et al., in prep.) corresponds to AzTEC42 in 
our nomenclature (cf. lScott et al. 20081 Table 1 therein). 
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Fig. A.l. continued. 


AzTEC17. The 1.3 mm source AzTEC17a (lying 473 
SW of the pointing centre, and detected at a significance 
level of 6.2<x) is clearly associated with the 20 cm source 
COSMOSVLADP_J095939.19+023403.6 (S 2 ocm = 68 ± 

13 juJy), and a VLA 10 cm source (37 pJy beam -1 or 8.2<x). 
The source also shows Spitzer /MIPS and IRAC emission. There 
is a Herschel 250 pm source 3705 east of AzTEC17a [ID 1753 
in the COSMOS SPIRE 250 pm Photometry Catalogue from 
HerMES (IQliver et al. 20121) 1. The search for Herschel counter¬ 
parts was performed by using a search radius of 971, i.e. half 


the SPIRE beam FWHM at A = 250 pm. AzTEC 17a has an 
optical-NIR counterpart about 174 SW of the 1.3 mm emis¬ 
sion peak [ID 1475165 in the COSMOS photometry catalogue 
( |Capak et al. 2007 1 )] . We note that the source visible in the Ultra¬ 
VISTA and ACS I -band images, lying 176 NW in projection 
from AzTEC 17a, is the galaxy COSMOS J095939.12+023405.1 
( |Capak et al. 2007 1 ), which has a photometric redshift of z = 
0.793 (the source 271694 in the DR 1 UltraVISTA photometric- 
redshift catalogue: lllbert et al. 20131) . Cross-correlation with the 
COSMOS photometry catalogue yielded a candidate optical 
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Fig. A.l. continued. 


counterpart for AzTEC17b (ID 1475223), about 172 SW of the 
PdBI emission peak. However, AzTEC 17b has no counterparts 
at UltraVISTA bands or at MIR or cm wavelengths. There are 
two negative sources (-4.5<x and -6.2<x) within < 1" of the 
phase centre. As AzTEC 17a is a confirmed SMG, this could 
mean that the 4.5<x source AzTEC 17b is spurious. 

For AzTEC17a, the primary photo-z solution is z = 0.75^ 
Because the source of this photo-z lies about ~ 174 from 
the PdBI position, it is questionable whether it is related to 
AzTEC 17a (although it is within the synthesised beam). The 
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primary photo-z value is however comparable to the very se¬ 
cure (quality flag 3) spectroscopic redshift z spe c = 0.834 mea¬ 
sured towards AzTEC 17a (0726 offset from the PdBI position) 
with Keck/DEIMOS (M. Salvato et al., in prep.). We note that 
for the measurements the slit was centred on a position 1760 
away from AzTEC 17a (see Appendix E), and that there are no 
HST /ACS I -band sources within the slit boundaries. For com¬ 
parison, the angular offset between the 1.4 GHz radio source 
and the PdBI detection peak is only 0722, and the redshift de¬ 
rived from the radio/submm flux density ratio is z = 2.29 + 0.42. 
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Fig. A.l. continued. 


A comparable redshift of z - 2.77 was derived from the Ward- 
low et al. [2011; their Eq. (1)] redshift estimator based on the 
Spitzer /IRAC 3.6 yum and 8 /im flux densities. The last two val¬ 
ues agree with a shallow “bowl” in the^ ot distribution at z - 2.7, 
but the z sp ec value is adopted in the present work. For AzTEC 17b 
we derive a photo-z of 4.14^ There is a dip in the^ ot distribu¬ 
tion also at z - 0.4. However, the 1.4 GHz non-detection towards 
AzTEC 17b results in a lower limit to its redshift of z > 2.49, 
which is consistent with the above photo-z value of 4.14^ ^. 

AzTEC18. The 1.3 mm 4.5cr source seen towards AzTEC 18 
has counterparts at optical and NIR wavelengths (0782 from the 


PdBI position), and is therefore considered a potential SMG. 
This is further supported by the fact that only one negative 
source, being of -4.2 cr significance, was detected in this field. 

AzTEC18 has a photo-z solution of z p hot = 3.00^ }^, which 
is consistent with the redshift derived from the radio/submm 
flux-density ratio of z > 2.20. 

AzTEC 19. AzTEC 19a lies 371 NE of the SCUBA-2 
source SMMJ100028.6+023201 (or 450.00 or 850.07) identi¬ 
fied by Casey et al. (2013). This angular offset is within the 
JCMT/SCUBA-2 beam size (FWHM) of ~ 7" at 450 /im. With 
the deboosted flux densities of S' 450 ^ = 37.54 ± 6.58 mJy 
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Fig. A.l. continued. 


and S 850 yum = 9.21 ± 1.45 mJy, this was the strongest 450 pm 
source found by Casey et al. (2013) in the COSMOS field. 
These authors identified two possible optical counterparts to 
SMMJ100028.6+023201 (see their Table 6), but our higher- 
resolution observations show that only one of them - ly¬ 
ing 1703 west of the PdBI peak position - can be taken 
as a candidate counterpart (the other source lies 2765 SW 
of the PdBI peak). AzTEC19a was also detected by Her- 
schel. In the COSMOS SPIRE 250 pm Photometry Cata¬ 
logue, the source ID is 2277 (0743 offset). AzTEC 19a is as¬ 
sociated with both a 20 and 10 cm radio-continuum source. 


In the VLA Deep Catalogue, the corresponding source has 
the name COSMOSVLADP_J100028.70+023203.7 (S 2 0cm = 
78 ± 12 juJy). The 10 cm peak flux density is Siocm = 
28.8 pJy beam -1 , making it a 6.4<x detection. The source is 
also associated with Spitzer IR emission. The 9.7<x 1.3 mm 
source AzTEC 19b lying at the border of the PdBI PB has some 
Spitzer /IRAC 3.6 pm and lower wavelength emission just north 
of it. This emission can be associated with the galaxy COSMOS 
J100029.24+023211.5 ( [Capak et al. 2007| , which has a photo-z 
of about 1.27 (source 262768 Gilbert et al. 20131) . As it lies 1767 
NW of the 1.3 mm peak, it is probably unrelated to AzTEC 19b. 
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Fig. A.l. continued. 


However, there is also a NIR source within the 3<x contour of 

1.3 mm emission, about 075 from the mm peak. In the COSMOS 
ACS /-band photometry catalogue (BLeauthaud et al. 2007]) . the 
ID of this source is 1486, while in the DR1 UltraVISTA cata¬ 
logue its ID is 262766, and its reported photo-z value is about 
1.30 (I lbert et al. 2013 ). Three negative sources (-4.8cr, -6.6<x, 
-9.5 cr) were detected in the AzTEC 19 field. However, only one 
of them (-6.6 cr) lies within the PB FWHM (979 from the phase 
centre), while the remaining two lie outside the PB (13" - 1378 
away from the phase centre). Because AzTEC 19a is confirmed, 
and AzTEC 19b has a high 1.3 mm detection S/N of 9.7 and is 
associated with multiwavelength emission, we are not expecting 
to have any spurious sources in this field. 

The optical/IR counterpart of AzTEC 19a (CFHT i* - K = 
1.49), located only 072 from the PdBI position, has a photo- 
z solution of z = 3.20^^. There is a spec-z value of 1.048 
measured for a source only 07195 from the PdBI position of 
AzTEC 19a with the VLT Visible Multi-Object Spectrograph 
(VIMOS) in the zCOSMOS project. However, the correspond¬ 
ing quality flag is 1.1, meaning that the z spe c value is insecure 
(< 25% reliability). For the 1.4 GHz source - situated 0753 from 
the PdBI position - we derived a radio/submm-based redshift 
of 4.22 ± 0.91, which is comparable within the errors with our 
photo-z solution. For comparison, Casey et al. (2013) derived 
a photometric redshift of z = 2.86^ ^ for their SMG source 
SMMJ100028.6+023201, which is associated with AzTEC 19a. 
A comparable value of z ~ 2.3 can be derived from the Spitzer 
IRAC/MIPS flux densities [Pope et al. 2006} their Eq. (2)], and 
the value z - 2.90 is obtained when using the IRAC 3.6 pm and 
8 pm flux densities (Wardl ow et al. 20 lTt . In the present paper, 
we adopt our photo-z solution of z = 3.20+® 4* for AzTEC 19a. 
For AzTEC19b, the photo-z solution of z P hot = 1.11 ± 0.10 is 
adopted as the redshift of the aforementioned SMG, while the 

1.4 GHz non-detection suggests a very high radio/submm-based 
redshift of z > 6.57, where no SMGs have been discovered to 
date. There is, however, a degeneracy between the dust temper¬ 


ature of the source and its redshift (both affecting the source 
SED), and the radio dimness of AzTEC 19b could in principle 
be the result of a low dust temperature (e.g. IB lain et aT~2 002 [ 
IKovacs et al. 20061) . 

AzTEC20. Interestingly, the 6 cr source about 578 north of 
the phase centre shows no emission at other wavelengths. The 
~ 3 cr features of 1.3 mm emission seen near the phase cen¬ 
tre appear instead to be associated with a source seen at sev¬ 
eral different wavelengths from 20 cm to NIR. The Ultra¬ 
VISTA catalogue ID of the latter source is 306331 (z - 1.98; 
I lbert et al. 2013 ). The 20 cm source has a peak flux density 
of 42.3 ply beam -1 , hence has a S/N ratio of about 3.5. We 
note that the VLA Deep catalogue contains sources down to 
4 cr or about 48 /Jy beam -1 ( Schinnerer et al. 2010 ). About 577 
west of this source, a slightly stronger 20 cm source candi¬ 
date (45 pJy beam -1 or ~ 3.8<x) can be seen. The 10 cm 
source near the PdBI phase-tracking centre is a ~ 6.2<x detec¬ 
tion (28.1 piy beam -1 ). The Spitzer /MIPS 24 pm emission is 
quite extended, but is clearly resolved into two sources in the 
IRAC 3.6 pm image and yet more sources in the UltraVISTA 
NIR im ages. The Herschel/ HerMES/SPIRE 250 pm catalogue 
dOliver et al. 20121) contains a source (ID 3076) near the phase 
centre, 772 from our PdBI source. There is one negative source 
of -4.8<t significance in this field, located on the SE side of our 
1.3 mm source. Although the positive source is more significant 
(6cr) compared to the negative feature, and it fulfils our detection 
criterion of S/N > 4.5, it has no optical-to-IR counterparts that 
could confidently confirm that it is real. The lack of radio emis¬ 
sion from AzTEC20 yields a lower redshift limit of z > 2.35. 

AzTEC21. The south-western clump of the detected filamen¬ 
tary structure, AzTEC21a, is associated with a 20 cm source of 
peak flux density of 63 pJy beam -1 (~ 3.9<x). This also coincides 
with the position of a Spitzer IR source. AzTEC21b is proba¬ 
bly part of the same structure (see below). The bright galaxy ly¬ 
ing 174 SE of AzTEC21c is COS MOS J100002.93+024639.9 
(V = 20.681; [Capak et al. 2007[ ). The ID of this galaxy in 
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the DR1 UltraVISTA catalogue is 327783, and its photo-z is 
about 0.34 ( lllbert et al. 20131) . We note that the source ID in 
the zCOSMOS catalogue is 846495, but its spectroscopic red- 
shift measurement could not be attempted (confidence class 0; 
[Lilly et al. 2007] 2009). 

For AzTEC21a, the photo-z solution is z = 2.60^ }^ (for a 
source with the CFHT colour i* - K = 2.39). The optical/NIR 
source lies only 0713 away (NE) from the PdBI source according 
to the previous COSMOS/UltraVISTA catalogue (ID 1688587), 
but in the most recent UltraVISTA-TERAPIX DR the nearest 
source (ID 328878) lies 171 NE, i.e. 8.5 times further away from 
our source, making the reliability of the proposed counterpart 
questionable. The 1.4 GHz non-detection (only 3.9<x) suggests 
a high redshift of z > 3.45. On the other hand, the Spitzer pho¬ 
tometric redshift is estimated to be about 1.5 ( [Pope et al. 2006| ), 
and the value z - 2 is derived using the Wardlow et al. (2011) 
IRAC flux-density method. The last value is comparable to our 
photo-z of z = 2.60^ }^. For AzTEC21b, we derived a red¬ 
shift of zphot = 2.80^-jg (source 0745 from the PdBI peak; 
CFHT//* - K = 1.74), while the radio/submm based value is 
z > 2.47, consistent with our photo-z solution. The most recent 
UltraVISTA DR does not contain a nearby (< 274) counterpart 
to AzTEC21b. We note that the overall 1.3 mm emission mor¬ 
phology (Fig. |A.l| ) could indicate a relation between AzTEC21a 
and 21b, and this is further supported by their comparable photo- 
z values of z = 2.60 +q and z = 2.80^ }^ (i.e. their redshifts 
could be identical). For AzTEC21c, we derived a lower limit of 
z > 1.93 from the upper limit to the 1.4 GHz flux density. 

AzTEC22. The candidate point-like 5.1<x PdBI source at 
the southern edge of the field has no counterparts at other 
wavelengths. The detection of two negative features with 
|S/N| = 5.7 and 5.9 in this field provides a hint that the 
1.3 mm detection is spurious. A submm source was de¬ 
tected with Herschel towards our phase centre (250 pm ID 
5470), about 8" north from the above mentioned PdBI fea¬ 
ture. Moreover, about 275 from our pointing centre, there is 
the 20 cm source COSMOSVLADPP_J095950.57+022827.5 
(*S 20 cm = 124+12 ply), which is also seen in the VLA 
10 cm image (58.6 ply beam -1 or ~ 13<x). The Spitzer IR 
images show a “double source”, and a stronger 24 pm emit¬ 
ter is associated with the above-mentioned radio-continuum 
source J095950.57+022827.5. A trace of 20 cm emission 
(43.8 ply beam -1 peak surface brightness) can also be seen to¬ 
wards the position of the weaker 24 pm source. The two Spitzer 
sources can be seen in the UltraVISTA NIR images: the NW one 
has the ID 244762 (Zphot - 1 .8 - 1.9 ), while the SE source is 
244405 at z p hot - 2 (lllbert et al. 20131) . For our PdBI detection 
the lack of a radio counterpart suggests a redshift of z > 3.0. 

AzTEC23. The 1.3 mm feature seen towards this source (574 
NW of the AzTEC centroid) has counterparts at optical-NIR 
wavelengths (1724 from the PdBI position). Cross-correlation 
with the Herschell SPIRE 250 pm catalogue shows the presence 
of a source (ID 2659) about 3" NW from the 1.3 mm feature. 
A visual inspection of the VLA 20 cm image reveals an EW- 
oriented, elongated emission feature, whose western emission 
peak (2748 from the 1.3 mm source) has a peak surface bright¬ 
ness of 43.5 ply beam -1 and the eastern peak has the surface 
brightness of 39.8 ply beam -1 . 

We derived a primary photo-z of z = 1.60^ ^ with a sec¬ 
ondary solution at z - 4.3, while the radio non-detection implies 
the lower limit z > 2.06. We adopt the redshift z = 1.60^ ^ but 
note that the Xtot distribution of AzTEC23 is quite complex. 
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AzTEC24 . This source is called AzTEC/C48 in the 
ASTE/AzTEC 1.1 mm catalogue of Aretxaga et al. (2011). The 
ASTE/AzTEC peak position lies 878 NE of the JCMT/AzTEC 
centroid. We have found three candidate PdBI sources of 4.9 - 
5.1<x significance (outside/at the border of the PB), but none 
of them have counterparts at other wavelengths. Altogether 
five negative sources with |S/N| = 4.3 - 5.9 were detected 
in the field, and two of them lie outside the PB FWHM. 
Hence, the identified positive 1.3 mm sources might be spuri¬ 
ous and should be treated with caution. About 575 SW from 
AzTEC24b there is a Herschel submm source (250 pm ID 4991). 
Moreover, 275 SW from AzTEC24b there is the VLA 20 cm 
source COSMOSVLADP_J100039.28+023845.3 (S 2 ocm = 63 + 
13 ply). Aretxaga et al. (2011) associated this source with the 
radio counterpart of their source AzTEC/C48 (277 separation). 
In the VLA 10 cm image, a source with a peak flux density of 
28.5 ply beam -1 (6.3<x) can be seen at the 20 cm source po¬ 
sition. The cm radio-continuum source is also associated with 
Spitzer IR emission. When referring to the COSMOS catalogue 
of Capak et al. (2007), the source can be identified as the galaxy 
COSMOS J100039.29+023845.4, and in the UltraVISTA cata¬ 
logue the source ID is 293896 (z p hot - 2.1: lllbert et al. 20l3l) . 
For the three radio non-detected components AzTEC24a, 24b, 
and 24c we derived lower redshift limits ofz > 2.35, z > 2.28, 
z > 3.17, respectively. 

AzTEC25. In this field, none of the 1.3 mm point-like fea¬ 
tures fulfilled our detection criteria. We note that 1737 o f the 
4.2 cr feature (the most south-western feature shown in Fig. |A. l[ ) 
there is the source 1925434 from the COSMOS+UltraVISTA 
catalogue, but its CFHT//*-band magnitude (24.27 + 0.49 mag) 
with respect to the CFHT/X-band and UltraVISTA X s -band 
magnitudes (23.56 + 0.20 mag and 24.30 + 0.23 mag, respec¬ 
tively) suggests a very blue colour index, hence it is unlikely re¬ 
lated to the 1.3 mm feature]^] There are also two negative sources 
(-4.3<x and 5.0<x) in the field, inside the PB FWHM. In the VLA 
20 cm map, a source candidate with a peak surface brightness of 
58.7 ply beam -1 (only ~ 3.3<x) is visible close to our pointing 
centre (1" away). However, no shorter-wavelength emission can 
be seen towards this source. 

AzTEC26. Two 1.3 mm point sources were detected around 
the phase centre of this source, and the source we called 
AzTEC26a has an optical-to-NIR counterpart (0794 from the 
1.3 mm peak). The weaker, 4.8<x source AzTEC26b has no mul¬ 
tiwavelength counterparts. The VLA 20 cm image detects some 
emission near both sources, but the S/N ratio of the radio detec¬ 
tion is less than 2. 

For AzTEC26a, we derived a photo-z value of z = 2.50+g^. 
Although the Subaru// + -CFHT/X colour of this source, 1.06, is 
relatively blue, the radio/submm-based redshift of z > 1.87 sup¬ 
ports the aforementioned z p hot solution. The radio non-detection 
of AzTEC26b gives a lower redshift limit of z > 1.79. 

AzTEC27. There is a hint of 20 cm emission associ¬ 
ated with the elongated 1.3 mm source, although its 20 cm 
peak surface brightness is only 32.1 ply beam -1 or about 
2.5 cr. No other wavelength counterparts are detected. At the 
projected distance of 572 NE of the 1.3 mm peak po¬ 
sition, there is the bright (V = 20.359) galaxy COS¬ 

MOS J100039.47+024055.5 ( [Capak et al. 2007| . The Ultra- 


9 If the dust column density in the SMG is high, it is expected to appear 
reddened. However, if the spatial distribution of dust is inhomogeneous, 
the unobscured star formation could cause the galaxy to appear blue. 
Moreover, the viewing angle can contribute to the observed effective 
colour of the galaxy (see e.g. |Wang et al. 2013] ). 























Miettinen et al.: Plateau de Bure Interferometer 1.3 mm imaging of SMGs 


VISTA catalogue gives a photo-z of about 0.25 for this 
galaxy (ID 303584; lllbert et al. 2001 About 2'.'9 SW of the 
1.3 mm peak position, there is another galaxy, namely COS- 
MOS J100039.07+024050.2 i ]Capak et aI. 2007) . whose photo- 
metric redshi ft is about unity (source 303782 in the UltraVISTA 
catalogue; lllbert et al. 2013i . As described in Appendix D, 
AzTEC27 appears to be subject to gravitational lensing by the 
foreground galaxies J100039.47 and J100039.07, and our lens 
model suggests that their combined lens effect amplifies the 
/fobs = 1.3 mm flux density by a factor of ~ 2. At the distance of 
475 north from AzTEC27, there is a Herschel/ HerMES 250 p m 
source (ID 1500: [Oliver et al. 2012 ). Based on its radio dimness, 
we derived a redshift of z > 4.17 for AzTEC27, which makes this 
source potentially the highest-redshift SMG among AzTEClb- 
30. 


AzTEC28. This clearly detected 1.3 mm source, 373 SE of 
the phase centre, has no counterparts at other wavelengths shown 
in Fig. A.l| Casey et al. (2013) detected AzTEC28 with SCUBA- 
2 (their source SMMJ100004.5+023042 or 450.20), with the 
450 p m peak lying about 572 NW of the PdBI peak posi¬ 
tion. The deboosted flux density at 450 /dm was reported to be 
19.11 ± 5.91 mJy. The optical counterpart (at z p hot = 0.76+gj^) 
reported by Casey et al. (2013) lies 4767 north of the PdBI 1.3 
mm peak, hence is unrelated to the SMG. We also note that the 
ASTE-detected source AzTEC/C 150 from Aretxaga et al. (2011) 
lies 8744 NE of our PdBI 1.3 mm source - still within the 34" 
beam of ASTE/AzTEC at 1.1 mm. About 676 NE of AzTEC28, 
there is the Herschel 250 pm source 4388 from the HerMES sur¬ 
vey dOliver et al. 2012t . 


AzTEC29. There are no clear signatures of PdBI 1.3 mm 
emission inside the PB. The source candidate AzTEC29a (4.7<x) 
lies at the border of the PB and the strong source AzTEC29b 
(7.3a-) at the northern edge of the map, is outside the PB. 
The latter could be associated (0776 offset) with the source 
1685295 from the COSMOS+UltraVISTA catalogue. Two neg¬ 
ative sources of -4.3<x and -5.4<x were detected in the field out¬ 
side the PB FWHM. Hence, AzTEC29a, which does not show 
up at other wavelengths, could be spurious. The VLA 20 cm im¬ 
age of the source region detects emission in its south-western 
corner (peak surface brightness of 50.5 pJy beam -1 or ~ 3.4<x). 
The Spitzer/ MIPS 24 jdm emission near the field centre can 
be associated with the galaxy COSMOS J100026.79+023749.4 
( [Capak et al. 2007|) or_J>89240 in the UltraVISTA catalogue 
(Zphot - 0.58: !Ilbert et al. 20131) . 

For AzTEC29a, we derived a radio/submm-based redshift 
of z > 2.96. The photo-z of 1.45+Q3g derived for AzTEC29b 
(Subaru// + -CFHT/AT = 2.67) is much lower than the unrealisti¬ 
cally high value z > 7.25 derived from the radio/submm flux- 
density ratio, hence the photo-z value is adopted. 


AzTEC30. A 4.6<x candidate point source with no counter¬ 
parts is detected in this field. There is a Spitzer/ MIPS 24 pm 
source in the NW part of the field, ~ 475 SE of our tentative 
source, that is also visible at the IRAC wavelengths but has no 
catalogue identification. The radio non-detection of AzTEC30 
gives a redshift of z > 2.51. 


Finally, we note that none of our sources are detected in 
X-rays, implying that none host a prominent AGN. Cross¬ 
correlation with the Chandra- COSMOS Bright Source Cata¬ 
logue v2.1 and COSMOS XMM Point-like Source Catalogue 
v2.0 revealed that the nearest X-source to any of our SMGs is 
XMMU J100002.8+024635, lying 9796 south of the PdBI phase 
centre towards AzTEC21. However, it is possible that some of 
the studied SMGs host an extremely Compton-thick AGN (i.e. 


with obscuration due to high column densities of dust) that re¬ 
mains undetected in the existing X-ray images. 

Appendix D: Gravitational lens modelling of 
AzTEC27 

There are two foreground galaxies seen in projection close to 
the SMG AzTEC27 (see Fig. |D.1| ). The north-eastern galaxy, 
at (Aar, AJ) = (+4702, +3736) from the 1.3 mm peak position 
of AzTEC27, is COSMOS J100039.47+024055.5. The photo¬ 
metric redshift and stellar mass ofJl 00039.47 are z p hot - 0-25 
and log(M*/M©) = 10.084 ( lllbert et al. 20131) . On the south¬ 
western side, at (Aar, Ac)) = (-1798, -1794) from AzTEC27, 
the foreground galaxy is COSMOS J100039.07+024050.2 at 
Zphot = 0.998 with a stellar mass of log(M*/M 0 ) = 10.713. 
AzTEC27 is potentially subject to gravitational lensing by these 
two intervening galaxies, and therefore to better understand its 
intrinsic physical properties requires a lens model. 


o 


J100039.47+024055.5 (z=0.25) 
0 



J100039.07+024050.2 (z=0.998) 


Fig. D.l. HST/ACS I -band (/Central = 8 333 A) image towards AzTEC27 
(linear colour scale), overlaid with white contours showing the PdBI 
1.3 mm emission. The positive contour levels are as in Fig. [l] The green 
dashed circles of radius 5724 (20.31 physical kpc at z = 0.25) and 2777 
(22.26 physical kpc at z = 0.998) are centred on the galaxies COS¬ 
MOS J100039.47+024055.5 and J100039.07+024050.2, respectively. 
The synthesised PdBI beam (1776 x 1702, PA. 5?19) is shown in the 
bottom left corner. 


To estimate the strength of the lensing effect for AzTEC27, 
we carried out a gravitational lens modelling using the publicly 
available Python software called uvmcmcfit[^]that will be de¬ 
scribed in detail by R. S. Bussmann et al. (in prep.). The code is 
a modified version of the one used in the papers by Bussmann et 
al. (2012, 2013) and uses the visibilities to determine the good¬ 
ness of fit. To sample the posterior probability density function 
of our model parameters, we used the MCMC sampling code 
emcee ( |Foreman-Mackey et al. 2013| . 

The lensed background source in the model is assumed to be 
an elliptical Gaussian source, described by the following param¬ 
eters: the total intrinsic (unlensed) flux density (5 m), the effec¬ 


tive radius (R e ff = Va x b, where a and h are the semi-major and 
semi-minor axes), the projected angular offset from the model 
image centroid, the axial ratio (b/a), and the position angle (P.A. 
measured E of N). The lens is assumed to be a singular isother¬ 
mal ellipsoid (SIE), parameterised by the angular Einstein radius 
(#e), the angular offset from the model image centroid, the axial 
ratio, and the P.A. The magnification factor is then computed as 


10 https://github.com/sbussmann/uvmcmcfit 
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p = S out /^ in j where Scut is the source’s total lensed flux density 
in the best-fit model. 

The modelling was performed assuming three different sce¬ 
narios: i) J100039.07 at z p hot = 0.998 is acting as a lens; ii) 

J100039.47 at z p h 0 t - 0.25 is responsible for lensing; and iii) 
both the above intervening galaxies act as lenses. The magnifi¬ 
cation factor in these three cases was found to be /i = 1.36+0.11, 

1.17 ± 0.04, and 2.04 ± 0.16, respectively. Both the galaxies 
J100039.07 and J100039.47 are therefore causing only a weak 
lensing effect, the former one, having a ~ 4.3 times higher stel¬ 
lar mass and lying at a higher redshift (i.e. closer to AzTEC27) 
than the latter, being the slightly more stronger lens. The weak 
lensing is consistent with the fact that we see only a single image 
of AzTEC27 in the PdBI 1.3 mm map. In the present paper we 
assume the two-lens system and adopt the value p = 2.04 ± 0.16 
for AzTEC 27. 


Appendix E: The DEIMOS spectrograph slit 
parameters 

In Table EH we list the central coordinates, sizes (length 
and width), and position angles of the slits used for the 
Keck/DEIMOS spectral line observations towards AzTEC2, 5, 
6, 9, 10, 13, and 17a. The DEIMOS slits are also illustrated in 
Fig. |E.1| where we show the Ultra VISTA 7-band NIR images 
towards the above sources. 
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Table E.l. Slit parameters of the Keck II/DEIMOS spectral line observations. 


Source 

<*2000.0 

[h:m:s] 

£ 2000.0 

Length 

["] 

Width 

["] 

P.A. 

[°] 

Offset" 

["] 

AzTEC2 

10 00 07.98 

+02 26 12.2 

7.7 

1.0 

90 

0.98 

AzTEC5 

10 00 19.69 

+02 32 04.4 

8.0 

1.0 

90 

0.90 

AzTEC6 

10 00 06.54 

+02 38 37.6 

7.2 

1.0 

90 

0.61 

AzTEC9 

09 59 57.28 

+02 27 30.6 

8.7 

1.0 

90 

0.45 

AzTEC 10 

09 59 30.88 

+02 40 33.9 

10.5 

1.0 

90 

1.80 

AzTEC 13 

09 59 37.01 

+02 33 20.0 

10.1 

1.0 

90 

0.60 

AzTEC 17a 

09 59 39.30 

+02 34 03.7 

9.4 

1.0 

90 

1.60 


Notes. (a) Angular offset between the SMG position and the slit centre (cf. Fig. E.l 



Fig. E.l. Ultra VISTA F-band (/t e ff = 1-02 /i m) images towards AzTEC2, 5, 6, 9, 10, 13, and 17a sh own in linear scale. All the images are 15" 
on a side. The white rectangles indicate the DEIMOS slit positions, sizes, and orientations (see Table |E.l| ). We note that all the slits were aligned 
horizontally along the east-west direction. The central position of the slit is marked with a white plus sign, while the green circle of radius 075 
shows the SMA 890 jam peak position f Younger et al. 2007} 2009) except in the case of AzTEC17a where it marks the PdBI 1.3 mm peak position. 
The red circle in the AzTEC6 panel represents the optical galaxy discussed by Koprowski et al. (2014); see Appendix B for details. 
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